Draft version November 20, 2012 

Preprint typeset using L^TgX style cmulateapj v. 5/2/11 



COSMOLOGY WITH PHOTOMETRICALLY-CLASSIFIED TYPE IA SUPERNOVAE FROM THE SDSS-II 

SUPERNOVA SURVEY 

Heather Campbell 1 , Chris B D 'Andrea 1 , Robert C. Nichol 1 , Masao Sako 2 , Mathew Smith 3,1 , Hubert Lampeitl 1 
Matthew D. Olmstead 8 , Bruce Bassett 4,5,6 , Rahul Biswas", Peter Brown 8,9 , David Cinabro 10 , Kyle S. Dawson 8 
Ben Dilday 11,12 , Ryan J. Foley 14,15 , Joshua A. Frieman 16,17,18 , Peter Garnavich 19 , Renee Hlozek 20 , 
Saurabh W. Jha 13 , Steve Kuhlmann 7 , Martin KuNz f1,21 , John Marriner 18 , Ramon Miquel 22,23 , 
Michael Richmond 24 , Adam Riess 25,26 , Donald P. Schneider 27,28 , Jesper Sollerman 29 , Matt Taylor 10 , 

Gong-Bo Zhao 1,30 

1 Institute of Cosmology and Gravitation, University of Portsmouth. Portsmouth, POl 3FX, UK 
2 Department of Physics and Astronomy, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104, USA 
3 Department of Physics, University of Western Cape, Bellville 7530, Cape Town, South Africa 
4 Mathematics Department, University of Cape Town, Rondebosch, Cape Town, South Africa 
5 South African Astronomical Observatory, Observatory, Cape Town, South Africa 
6 African Institute for Mathematical Sciences, Muizenberg, 7945, Cape Town, South Africa 
7 Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA 
8 Department of Physics & Astronomy, University of Utah, 115 South 1400 East 201, Salt Lake City, UT 84112, USA 
9 George P. and Cynthia Woods Mitchell Institute for Fundamental Physics & Astronomy, Texas A. & M. University, Department of 

Physics and Astronomy, 4242 TAMU, College Station, TX 77843, USA 
10 Wayne State University, Department of Physics and Astronomy, Detroit, MI 48126, USA 
11 Las Cumbrcs Observatory Global Telescope Network, 6740 Cortona Dr., Suite 102, Goleta, California 93117, USA 
12 Department of Physics, University of California, Santa Barbara, Broida Hall, Santa Barbara, California 93106-9530, USA 
13 Dept. of Physics and Astronomy, Rutgers University, 136 Frclinghuyscn Road, Piscataway NJ 08854, USA 
14 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 USA 

15 Clay Fellow. 

16 Department of Astronomy and Astrophysics, The University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 
17 Kavli Institute for Cosmological Physics, The University of Chicago, 5640 South Ellis Avenue Chicago, IL 60637, USA 
18 Center for Particle Astrophysics, Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, IL 60510, USA 
19 Department of Physics, University of Notre Dame, Notre Dame, IN 46556, USA 
20 Department of Astrophysical Sciences, Priceton University, Princeton, New Jersey 08544, USA 
21 Departemcnt dc Physique Theorique and Center for Astroparticle Physics, Universite de Geneve, Quai E. Anscrmet 24, CH-1211 

Geneve 4, Switzerland 

22 Institut de Fi'sica dAltes Energies, Universitat Autonoma dc Barcelona, E-08193 Bellaterra (Barcelona), Spain. 
23 Institucio Catalana de Recerca i Estudis Avangats, E-08010 Barcelona, Spain 
Physics Department, Rochester Institute of Technology, Rochester, NY 14623, USA 
25 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
26 Department of Physics and Astronomy, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, USA 
27 Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 
28 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 
29 The Oskar Klein Centre, Department of Astronomy, AlbaNova, Stockholm University, 106 91 Stockholm, Sweden and 
30 National Astronomy Observatories, Chinese Academy of Science, Beijing, 100012, P.R.China 

Draft version November 20, 2012 



ABSTRACT 

We present the cosmological analysis of 752 photometrically-classified Type la Supernovae (SNe la) 
obtained from the full Sloan Digital Sky Survey II (SDSS-II) Supernova (SN) Survey, supplemented 
with host-galaxy spectroscopy from the SDSS-III Baryon Oscillation Spect roscopic Survey (B OSS). 
Our photometric-classification method is based on the SN typing technique of lSako et al.l (|2011[ ). aided 
by host galaxy redshifts (0.05 < z < 0.55). SNANA simulations of our methodology estimate that 
we have a SN la typing efficiency of 70.8%, with only 3.9% contamination from core-collapse (non-la) 
SNe. We demonstrate that this level of contamination has no effect on our cosmological constraints. 
We quantify and correct for our selection effects (e.g., Malmquist bias) using simulations. When fitting 
to a flat ACDM cosmological model, we find that our photometric sample alone gives fl m = 0.24+q'q5 
(statistical errors only). If we relax the constraint on flatness, then our sample provides competitive 
joint statistical constraints on Q m and comparable to those derived from the spectroscopically- 
confirmed three- year Supernova Legacy Survey (SNLS3). Using only our data, the statistics-only 
result favors an accelerating universe at 99.96% confidence. Assuming a constant wCDM cosmological 
model, and combining with H , CMB and LRG data, we obtain w = — 0.96+o;}q, O to = 0.29lo!o2 an d 

flk = O.OO^q qj (statistical errors only), which is competitive with similar spectroscopically confirmed 
SNe la analyses. Overall this comparison is reassuring, considering the lower rcdshift leverage of the 
SDSS-II SN sample (z < 0.55) and the lack of spectroscopic confirmation used herein. These results 
demonstrate the potential of photometrically-classified SNe la samples in improving cosmological 
constraints. 

Subject headings: supernovae: general — cosmology: observations — surveys — distance scale 
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1. INTRODUCTION 

Supernovac (SNe) have historically been classified 
based on their optical spectroscopic properties (e.g., 
lFilip P enkolfT99l . Type la supcrnovae (SNe la) are dis- 
tinguished from other classes of SNe by their lack of Hy- 
drogen and Helium spectral features, and the presence of 
other spectral features such as Si II absorption at rest- 
frame wavelength 6150 A. This particular optical classifi- 
cation is unique, as it efficiently separates two distinct SN 
physical processes. The progenitors of SNe la are White 
Dwarfs (WDs), unlike other common categories of SNe, 
which result from the core collapse of stars with initial 
mass M > 8 Mq. As the progenitors for core-collapse 
(non-la) SNe, such as Type II SNe, are more luminous, 
the progenitor star has been identified from archival im- 
ages on multiple occasions (|Sma rtt 2009]), while there 
have not yet been any direct observations of the WD 
progenitor of SNe la (th ough deep limits do exist on 
SN 2011fe: llJeral]|20ll . 

However, the lack of H and He in the spectrum, 
the composition of the ejecta, and the energy released 
in the explosion all strongly indicate that SNe la are 
the visible manifestation of a thermonuclear runaway ex- 
plosion in a carbon-oxyg en WD as its mass approaches 
the Chandrasekhar limit dHoyle fc Fowlerlll960l : iNomotol 
[1981 llben fc Tutukovl I1984D . The exact nature of 
the binary progenitor system (a single degenerate ob- 
ject accreting mass from a companion, or the merger 
of two WDs) has long been an open question. Re- 
cent observation s have shown that both channe l s can 
lead to SNe la dNueent et all 1201 fl iBloom et all 120111 : 
iDildav et al.l I2012D . making the relative prevalence of 
each channel of primary concern amongst progenitor 
studies. 

For over two decades SNe la have been of inter- 
est to cosmologists as distance indicators, as they dis- 
play less dispersion at their peak magnitude than other 
classes of supernovae and have high optical luminosities 
(> 10 9 L ). However, SNe la are not simple "standard 
candles", and the usefulness of SNe la are greatly en- 
hanced by our ability to standardize their magnitudes. 
The discovery that photometr ic properties of SNe la, 
such as the light-curve width (jPhillipsI Il993| ) and color 
(|Riess et al.lfl996l PTvTBpI Il997l ) . arc correlated with the 
absolute magnitude at peak allowed for accurate dis- 
tance measurements using SNe la, reducing the disper- 
sion in the measured distance modulus (fi) to ~0.14 mag. 
This work led directly to the discovery that high-redshift 
SNe la appear fainter than expected unless the expan- 
sion of the universe i s accelerating (|Riess et al.l 119981 : 
iPerlmutter et al.lll999D . 

Over the last decade SN surveys have evolved, span- 
ning greater fractions of the sky and discovering SNe 
at higher redshifts, with better photometric calibra- 
tion and an improved understanding of systematic un- 
certainties. Due to the differing observational require- 
ments necessary for monitoring SNe la over a wide 
range of redshifts (survey depth, area, wavelength cov- 
erage, etc.), the redshift-distance relationship, or "Hub- 
ble diagram", for SNe la is comprised of data from a 
number of surveys. The Hub ble Space Telesco p e Prq - 
gram (GOODS SN sample; IRiess et al.1 l200i 120071) . 
Supernovae Cosmology Project (SCP or HST Clus- 



ter SN sample: iKnop et al.l 120031: iKowalski et al.l 120081: 
lAmanullah et al.l 120101: iSuzuki et al.l 120111) the Super- 
nova Legacy Survey (SNLS: lGuv et al.ll2010f) , and Equa- 
tion of Stat e: SupErNovae trace Cosmic Expansion 
(ESSENCE; iMiknaitis etUl [2001 provide nearly all 
SNe la measurements at z > 0.4, with the first two sur- 
veys dominating at z > 0.9. 

The Sloan Digit al Sky Survey II Supernova Survey 
(SDSS -II SN survey; iFrieman et all[200a IKessler et all 
l2009aj ) has populated the Hubble diagram at interme- 
diate redshifts (0.1 < z < 0.4). Panoramic Survey 
Telescope & Rapid Response System (PanSTARRS) is 
currently finding thousands of SNe la in the same red- 
shift range. The largest of the low redshift surveys in- 
clude the Harvar d-Smithsonian Cent er for Astrophysics 
SN group (CfA; iHicken et a , J I2009I). 



the 



Nearby Su- 
pernova Factory (SNfactory; lAldering et al.1 l2002n. the 
Carnegie Supernova Project (CSP: lHamuv et alTl2006l : 
iFolatelli et al.l I2010D. the Lick Ob servatory Supernova 



Search ( Ganeshalingam et al.ll2010D . the Palomar Tran- 



sient Factory ([Law et a-1 .1120 091)7 and the C atalina Real- 
Time Transient Survey fCRTS: lDrake et al.l[2009[ ). Using 
a combination of data from different SN surveys, span- 
ning the full range of redshifts available, the most re- 
cent cosmological analysis using only SNe la finds tt m — 



0.18l£^ (stats) ±0.06(sys) and w = -OJOrJ^stats) 
^ail( s y s ) ([Conlev et al.ll20lH ). under the assumption of 
a flat Universe and a constant equation-of-state of dark 
energy (u>). 

A common theme across most previous SN cos- 
mology surveys is that SN spectroscopy is used only 
to classify the SN event and obtain a redshift (with 
the exceptio n of the SNfactory and, in some cases , 
the SNLS (IBronder et al.l [20071: IWalker et all I20T0I: 
lEllis et all 120081 )). The distance to each object clas- 
sified as a SN la is then determined from the multi- 
color time series photometry using o ne or more light - 
curve fittin g models, e .g., SA LT2 (IGuv et all 120071), 
MLCS2k2 (IRiess et al. | [1996 1 IJha et al.1 120071) . and 



SiFTO (|Conlev et al.l 120081) . Methods for constraining 
the absolute magnitude of an observed SN based on spec- 
troscopic line ratios have also been developed and ap- 
plied to nearby SNe (llblev et alJl2008t iBailev et al.H2009b 
Chotard et ail 120111: iNordin et al.1 120111 : iBlondin et al.l 



2012t iSilverman et al 120121 ). though as of yet they have 



not been tested over a cosmologically interesting redshift 
range. 

In the future, however, taking spectra of a large 
number of high-redshift SNe for classification pur- 
poses will be challenging, and could potentially limit 
the size and usefulness of any new survey. This 
will be the case for the Large Synoptic Surv ey Tele- 
scope (jLSST Science Collaboration et al.l I2009D . where 
the complete spectroscopic classification of all its SN can- 
didates will be simply impossible, thus other methods 
need to be employed to utilize these huge SN programs. 
This challenge is already being confronted by the latest 
SN surveys, such as the PanSTARRS, who have devised 
new and innovative techniques for classifying their SN 
candidates using only photometric imaging data (thu s 
far tested on a small subset; see IScolnic et al1[2009T) . 
The SN survey from the Dark Energy Survey (DES) will 
probably not be able to spectroscopically classify all of 
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their expected ~ 4000 high-quality, high-rcdshift SNc la 
(|Bernstein et alll2011l ). 

An alternative is to use photometr ic-only cla s sifica - 
tion techniques. This idea is not new - iPskovslril ()1977l ) 
proposed classifying SNe based on their observed de- 
cline rate - and even early SN la cosmology results in- 
cluded a significant fraction of high-r edshift even t s that 
lacked spectroscopic i d entification ( Riess et all 119981: 
Perlmutter et al.l 119991: iTonrv et al.l 120031: iRiess et al.l 



20041 ). One of the primary scientific drivers for de 



veloping photometric classifiers has been to aid in the 
spectroscopic follow-up of SN surveys, thus allowing 
them t o use their spectroscopic resources more effi- 
ciently (jSullivan et al.ll2006l ISako et al.ll2008D . However, 
making a Hubble diagram solely from photometrically- 
classified SNe requires a lower false-positive rate (i.e., 
contamination by non-la SNc) than does spectroscopic 
target selection from photometrically-classified SNe. 
Photometry-onl y Hub ble diagrams were introduced by 
iBarris fe Tonrvl (120041). and h a ve be en pre sented more 
recent ly by iRodnev fe Tonrvl (|2010f ) and iBazin et~aH 
(|20T1 . 

Most photometric classification methods fit ob- 
served light-curves to templates of different SN types 
and determine the likelihood o f each class. These 
methods ( e.g.. IPoznanski et all 120021: iSullivan et all 



2001 iJohnson fe CrottsT 120061: IPoznanski et al.l 



Kunz et all 



2007 



2007 



Kuznetsova fe Connolly 200., 

Rodney fe Tonrvl 120091: IGong et all 120101: iFalck et all 
20101: ISako et al.l 120111) typically remove SNe which 



resemble non-la SN templates based on their likelihoods, 
althoug h there is co n sidera ble variety in the details- 
Whil e iHlozek et all (|201lD (ba sed on iKunz et all 
20071 and further developed by INewling et all 120111: 
Knights et all I2012D uses templates to compute likeli- 
hoods for the type of each SN, they do not remove any 
objects from their cosmological fit. Rather, they use the 
likelihoods of each SN as a weight, retaining all possible 
information while computing an unbiased cosmology in a 
Bayesian manner. Higher level statistical analyses have 
also been applied to this problem, with a goal of finding 
a lower dimensional parameter space where there exists 
a cleaner separation between the different types, thus 
simplifying the classification problem. Examples of these 
approaches include se mi-supervised learnin g techniques 
such as diffusion maps ([Richards et al. 2012[), and kernel 
Principal C omponent Analysis (PC A) applied to SN 
light curves (llshida fe de Souzall20Tl . 

In anticipation of t he SN typing requir ements to 
be encountered by DES, iKessler et all (|2010aL hereafter 
KlOa) issued the "Photometric Supernova Classification 
Challenge" , providing simulated light-curves of different 
SN types based on a realistic DES-like SN survey and a 
training sample where the true SN type w as given. The 
results of this challenge are presented in IKessler et all 
(|2010bl hereafter KlOb), and provide some interesting 
insights into the relative performance of different SN clas- 
sifiers. Overall, several different classification strategics 
produce similarly high scores in terms of both efficiency 
and contamination, but all proved subject to significant 
level of contamination (~ 20%). A problem which is 
common for many methods is that they require a train- 
ing set of known SN types. If this set is biased and not 
representative of the whole sample then the classification 



will be biased as well; an effect that was seen in the KlOa 
challenge. 

In this paper we bu ild upon t h e ph otometric- 
classification algorithm of ISako et all (|201lL hereafter 
Sll), which obtained the highest overall Figure-of-Merit 
in KlOb. We use here the full three-year data-set from 
SDSS-II SN Survey, including a new collection of host- 
galaxy redshifts obtained by the SDSS -III Baryon Os- 
cillati on Spectroscopic Su rvey (BOSS; lEisenstein et all 
IMltlDawson et allboil ). 

We optimize selection cuts and determine the bi- 
ases of our new method with extens ive simulations using 
the SuperNova ANAlysis (SNANA; IKessler et all2009bD 
software package, and apply redshift-dependent correc- 
tions to our data. We show that photometric classifi- 
cation can provide SN la samples with low contamina- 
tion and well-understood biases, and present cosmolog- 
ical constraints that are competitive with those derived 
from existing spectroscopic samples. 

The paper is organized as follows. In Section [2] we 
detail the SN and host-galaxy data that we analyse in 
this paper. In Scction[3]we discuss our SN classifier, with 
emphasis on the light-curve fitter and selection criteria. 
We perform a rigorous analysis of and derive corrections 
for biases introduced by our selection criteria, shown in 
Sectional In Section [5| we present our full photometric 
Hubble diagram and consistency checks. Cosmological 
constraints from our photometrically derived sample are 
in Section [6] along with comparisons to other spectro- 
scopic cosmological fits. In Section [7] we discuss our re- 
sults, how this analysis could be improved, and how our 
work applies to upcoming large-scale SN surveys. Finally 
in Section[8]we detail the main conclusions of this paper. 

2. DATA 

2.1. The SDSS-II Supernova Survey 

The SDSS-II SN Survey is a dedicated search for 
intermediate-redshift SNe from repeated scans of the 
equatorial "Stri pe 82" region (cov e ring ~300 deg 2 ) of th e 
original SDSS (lYork et all f200Tl iFrieman et all 12008^ . 
For three months a year (Sep-Nov) ov er a three year pe- 
riod (2005-2007) the SDSS telescope (IGunn et alll2006l: 
IGunn et"all |1998|) perfor med multi-color uqriz imag- 
ing jFukugita et allll99l IGunn et alll2006Hlvezic et all 
120071 : iDoi et al.l l2010h ~of this area of sky, with a ca- 
dence of a few times per week. T he SDSS uses as- 
inh magnitudes (jLupton et al.l 11999^ . although SDSS is 
on the AB system after applying small offsets (see Ap- 
pendix 12.2. ip. The a nalys is of SDSS astrometry is de- 
scribed hilPkriFil] (|2003l) . 

This multi-epoch data was then used to iden- 
tify SN la candidates i n real-time for further 
spectroscopic observations (jSako et all 120081 ). result- 
ing in over 500 s pectroscopicall y-confirmed SNe la 
(iZheng et all 120081: lOstman et all 120101: iKonishi et all 
120111: iFolev et al.l I2012TT Well-observed subsamples 
of these intermediate-redshift SNe have been used in 
a variety of stu dies, primarily focus e d on constrain- 
ing c osmology (IKessler et all l2009ab iSollerman et all 
I2009t iLampeitl et al l l2010al). the mea surement of SN 
rates (jDildav et al.l 120101: iSmith et al.l 120111 ) and the 
study of host -galaxy properties an d their correlations 
with SNe la (jLampeitl et al.l l2010bl: iGupta et ail 120111: 



4 

ID 'Andrea etH] I20TH: iGalbanv et all 120121) . The full 
three-year SDSS-II SN Sample will be published in Sako 
et al. (in prep.). 

The real-time spectroscopic sample of SDSS-II SNc 
is incomplete and potentially biased, as a function of 
rcdshift. This bias comes from a number of different, 
and sometimes competing, effects and is therefore hard 
to predict a-priori and thus correct for. First, the deci- 
sions made by observers following-up SN candidates was 
based on the local weather conditions (at a variety of 
telescopes), the position of the SN candidates on the sky, 
and the location of the SN ca ndidate in the host galaxy. 
This is illustrated in Table 2 of lSmith et all (|201lD where 
the spectroscopic completeness of the SD SS-II SN Sur- 
vey dr ops below 40% at z > 0.4 (see also iKessler et all 
I2009al) . 

Second, our targeting of SN candidates gave prior- 
ity to_events i n red elliptical host galaxies (see Eqn 7 
of ISako et al.l (|2008f )). as these SNe are likely to be less 
affected by dust in their host galaxy. This prioritiza- 
tion is seen in the spectroscopically-confirmed SNe la, 
which have lower reddening values than predicted from 
simulations without spectroscopic selection (Figure 16, 
IKessler et al.l l2009aj ). However, the rate of SNe la 
in red elliptical galaxies is lower than seen in blue, 
star-forming spiral gala xies (|Mannucc i 2005; Wan g et all 
119971: iSmith et al.ll20lil) . so this additional up-weighting 
given to the SNe la in red ellipticals might be a subdom- 
inant effect. As will be shown in Figure [20]of this paper, 
the host galaxies of spectroscopically-confirmed SNe la 
are, on average, representative of the whole population of 
host -gal axy colors studied in our BOSS sample (see Sec- 
tion [472J . However, it is clear that the spectroscopically- 
confirmed SNe la from the SDSS-II SN Survey are bi- 
ased in absolute magnitude with respect to the entire 
SN la population. To avoid any such biases, we do not 
use spectroscopic SN information anywhere in our pho- 
tometric classification (as discussed in Section [3]). Any 
spectroscopically-confirmed SNc la that fail our photo- 
metric criteria are not included in our final sample to 
preserve selection consistency. 

2.2. BOSS Ancillary Targets 

We u ndertook a BOSS an cillary spectroscopic pro- 
gram (see iDawson et all (|2012l ) for details of these pro- 
grams), to obtain host galaxy information for most 
SDSS-II SN candidates. We aim to obtain a larger, and 
more complete, SN la sample from the SDSS-II SN Sur- 
vey through photometric classification. We also took 
the opportunity to target the host galaxies of a range of 
other, possibly interesting, transient events detected as 
part of the SDSS-II SN Survey. Such a project was well- 
suited to the small number of ancillary targets available 
in each BOSS spectroscopic field and would have been 
impossible to achieve on a normal instrument (either 
through queue-scheduling or normal observer-mode) , be- 
cause of the combination of the low surface density of 
faint targets (between 5 to 25 targets per deg 2 to r ~ 22) 
spread over a wide area (~ 3 00 deg 2 of "Stripe 82" ) . Ta- 
ble 1 of lFrieman et~aT1 IpOOl shows that the SDSS-II SN 
Survey used over 1000 hours (or ~ 100 nights) of tele- 
scope time between 2005 and 2006 to spectroscopically 
confirm over 350 SNe la, utilizing many of the larger op- 
tical telescopes in the world (HET, NTT, NOT, APO, 



Subaru, WHT, SALT, Keck, NOAO, TNG). It is hard to 
envisage how we could have used such resources to target 
thousands of host galaxies as presented herein. 

In detail, we obtained spectrosc opic observations us- 
ing the BOSS spectroscopic system l|Smee eTaIll20ll of 
3761 galaxies spread almost evenly across the "Stripe 82" 
region. These targets were chosen using two selection al- 
gorithms (described below) which were complementary, 
but different, in their scientific objectives. The first al- 
gorithm focused on improving cosmological constraints 
from SNe la in the SDSS-II SN Survey by obtaining a 
sample free from the possible spectroscopic-selection bi- 
ases discussed above. The other algorithm targeted in- 
teresting subsamples of transients detected as part of the 
SDSS-II SN Survey. In this paper we focus exclusively 
on the first of these objectives, but include in our analy- 
sis galaxies, and their associated SNe, observed by BOSS 
originating from cither algorithm. 

2.2.1. Algorithm One: Additional Type la candidates 

We targeted galaxies that hosted a SN event of any 
type, based on object classifications using the "Photo- 
metric SN IDentification" (PSNID)method of Sll, which 
we describe in detail in Appendix lAl We applied PSNID 
to the multicolor (ugriz) light-curve data from SDSS- 
II created with the Scene-Modeling Photometry (SMP) 
method of iHoltzman et ahl (|2008f P1. assuming a flat prior 
on all SN parameters. We do not include any spectro- 
scopic information for these objects, and thus place a 
flat prior on the SN redshift as well. We select transient 
events that were classified as likely SNc of any type, i.e., 
probable SN la, SN II or SN Ibc. 

For all these candidates, we visual ly inspected the 
SPSS images (from the DR7 Skyserver; lAbazaiian et al.1 
120091 ) of the three nearest galaxies to each transient and 
manually assigned the most likely host galaxy for each 
candidate. In the majority of cases (95%) the nearest 
galaxy on the sky (angular separation) to the candidate 
was classified as the host, but in some cases the near- 
est galaxy was either clearly a background object or a 
star misclassificd as a galaxy. In these cases we classi- 
fied either the second (4%) or third (1%) nearest galaxy 
as the most likely host. Taking into account the ob- 
servational limitations of the program, we gave priority 
to host galaxies with a fiber magnitude brighter than 
ffibcr = 21.25, based on SDSS-I/II photometry How- 
ever, to fully utilise our allocation of BOSS fibers, we 
also include a small subsample of host galaxies fainter 
than this limit. Combined, these samples made up our 
main target list of 2781 galaxies. 

In 66 cases, the nearest object to the SN event was 
classified as a faint star in the SDSS DR7 database, but 
was clearly a galaxy as seen in the SDSS images. For 
these objects, we targeted the nearest "stellar" source 
in the DR7 database, which could have a fainter fiber 
magnitude than our main BOSS targets (r^ber < 21.25) 
and were given lower priority for BOSS observations. 

1 Improvements in this photometry compared to that used in 
IKessler et al . (2009a) include the discovery of a small (~ 10%) cor- 
rection for an underestimate of uncertainties for low-flux data, and 
new (December 2011) AB offsets from the SDSS native magnitudes 
system of Au = -0.066, Ag = 0.021, Ar = 0.005, At = 0.020, and 
Az = 0.013. Details of this re-calibration can be found in Bctoulc 
et al. (in prep.) and Sako et al. (in prep.). 
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Finally, we cross-referenced the whole target list 
with the SDSS DR7 database and found 276 of these 
host galaxies had a spectrum already. For these cases we 
targeted the location of the SN event rather than placing 
the BOSS fiber at the centre of the galaxy, as was done in 
all other cases and in the original SDSS-I/II survey. The 
motivation for this fiber placement is that galaxy spec- 
tra at the location of the SN could be useful for studies 
into correlations between SN properties and the environ- 
ment in which the SN occurs ([Gallagher et al.ll2005l[200l 
Gupta et all 1201 li ID 'Andrea et all bo lit iGalbanv et al.l 
2012J). As our focus here is only on obtaining rcdshifts to 
aid in classifying transients, these spectra (at the location 
of the SN event) serve our purposes equally well. In these 
cases, we ignored the fiber magnitude of the host galaxy 
(calculated at the center of the galaxy) and observed the 
SN location position regardless of our r^bcr = 21.25 limit 
for the main sample. 

2.2.2. Algorithm Two: Random sample of additional 
transients 

The goals of our second sample were both to study 
our overall selection biases (e.g., determine the effect 
of AGNs) and to further study interesting and unusual 
variable objects observed by SDS S-II e.g., hydrogen- 
poor, superluminou s supernovae ([Quimbv et al.l 120111 : 
iLeloudas et al.ll2012[ ). Achieving the first of these goals 
requires an unbiased sample of non-SN transient host 
galaxies, created by choosing at random from the set of 
non-SN transients in the magnitude range of 19.5 < r < 
21.5. We imposed no magnitude limit on the host galax- 
ies of these targets, of which there were 980; however, 
we did require that the galaxy was detected in the DR7 
galaxy catalogue. 

As with the SN host-galaxy sample in Section [2. 2. 11 
we visually inspected the three nearest galaxies in DR7 
to each transient, selecting as our BOSS target the most 
likely host to the transient. This proved to be the nearest 
galaxy (angular separation) in the vast majority (99.4%) 
of cases. The high percentage of host galaxies being 
matched to the nearest galaxy, compared to the previous 
sample in Section 2.2.1, is likely due to these transients 
being quasars that are located in the cores of galaxies, 
as opposed to SNe, which can be located throughout the 
galaxy. The targets from this sample were given the low- 
est priority for observation when assigning BOSS fibers. 

Another key difference between Algorithms One and 
Two was that in Algorithm Two transient events were 
allowed to show variability over multiple years, or have 
light curves that failed in the initial photo-typing. Tar- 
gets selected by Algorithm One were detected in only one 
season of the SDSS-II SN Survey. 

2.3. Reduction of BOSS Spectra 

Our ancillary targets described above were merged 
with other ancillary BOSS targets on the "Stripe 82" 
region and observed together as part of the normal 
SDSS-III observing p r ogram during 2009 and 2010 (see 
lEisenstein et al.ll2011b iDawson et al.ll20T2T for a descrip- 
tion of BOSS observations and programs). Some of our 
targets were lost at this stage because of BOSS fiber colli- 
sion issues, i.e., two objects within 55 arc-seconds of each 
other cannot be observed on a single BOSS spectroscopic 
plate. 



By the end of 2010 all BOSS plates on "Stripe 82" , 
including most of our ancillary targets, had been ob- 
served. For the purpose of this paper we are only in- 
terested in the galaxy redshift measurements (and their 
errors). The details of the spectral reductions of our 
sample of BOSS SN host galaxy spectra can be found 
in Olmstead et al. (in prep.), and the redshift and ob- 
je ct classificatio n of BO SS spectra in general is described 
in iBolton et al.l (|2012f) . In brief, our sample was pro- 
cessed using the standard BOSS spectroscopic analysis 
software, which is based on the original SDSS-I/II reduc- 
tion pipelines. This pipeline has at least a 95% success 
rate in obtaining re dshifts for spectra fro m the primary 
galaxy sample (e.g.. lAnderson et al.ll2012l ). Olmstead et 
al. (in prep.) carried out a detailed comparison between 
the BOSS pipeline redshifts and those produced, and 
manually inspected, using the publicly available AAO 
RUNZ software. All the spectral data from our SN host 
galaxy ancillary progra m are now public as part of the 
SDSS DR9 data release (|Ahn et al.ll20T2h . 

After removing spectra with low redshift confidence, 
low signal-to-noise (S/N), and large redshift errors, our 
SDSS-II SN ancillary program on BOSS produced 3323 
reliable redshifts. This sample is composed of 2382 likely 
SN host galaxies from the 2781 targets selected via Al- 
gorithm One in Section |2.2.1[ and 941 from the random 
sample of 980 galaxies chosen in Scction [2.2.2l The spec- 
troscopic target efficiency from Algorithm One is 86%, 
which is approximately 10% lower than for Algorithm 
Two. This difference is due to two subsamples in Al- 
gorithm One that have a lower efficiency, probably due 
to the lack of an imposed fiber-magnitude limit in these 
cases. The first of these subsamples targets the near- 
est photometric object to the SNe when there was some 
ambiguity about the star/galaxy separation; this had an 
efficiency of only 27.3%, but was a very small subsam- 
ple. The second subsample was created using less strin- 
gent cuts on the quality of the light curve (i.e., lower 
S/N) to provide an additional list of probable SN loca- 
tions; this subsample has a lower average host galaxy 
fiber magnitude ((rfn, er ) = 22.17) than the main sam- 
ple ((rfiber) = 20.62). which leads to a lower redshift 
efficiency (70.1%). 

At this point, we checked the SDSS DR8 database 
for any host galaxies that failed to gain a redshift from 
our BOSS observations and reductions, based on either 
Algorithm One and Two as described above. We found 
an additional 178 host galaxies had a successful redshift 
measurement in DR8, and add them to our final BOSS 
sample. We note that most of these galaxies were on our 
BOSS target list, but were not observed during our an- 
cillary program. This combined sample (which, for sim- 
plicity, we will continue to call our 'BOSS' sample) forms 
the basis for our subsequent re-analysis of the SDSS-II 
SN light-curve data in Section [31 now with the SN can- 
didate redshift constrained to match the observed host 
galaxy. 

We present in Figure Q] the percentage of success- 
ful redshift measurements (from comparing with RUNZ 
and visual inspection) obtained from our BOSS-observed 
host galaxies, shown as a function of the host-galaxy 
r— band fiber magnitude. Figure [5] shows the measured 
redshift distribution for our BOSS targets, grouped by 
their target algorithm. For comparison we also show 
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Fig. 1. — The percentage of well-measured (based on comparison 
with RUNZ and visual inspection) redshifts obtained as a function of 
the ) — band fiber magnitude from the photometry of the SNc (Algo- 
rithm One) and general transients (Algorithm Two) host galaxies. 
We recover 3323 well-measured redshifts from the 3761 galaxies 
observed by our BOSS program. 
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BOSS Redshift 

Fig. 2. — Redshift distribution of the targets observed with 
BOSS. The black histogram is the subset which were selected 
using Algorithm One (candidate SN host galaxies) and the blue 
are the from Algorithm Two (host galaxies for general tran- 
sients) as discussed in Section [2] The green histogram shows all 
spectroscopically-confirmed SNe la from the SDSS-II SN Survey. 
The red histogram shows the additional 177 SDSS-II DR8 host 
galaxy redshifts. 

in Figure [2] the spectroscopically-confirmed SNe la from 
SDSS-II, which clearly peak at lower redshifts. There are 
a significant number of transient host galaxies with red- 
shifts greater than 0.5, extending out to z w 1, which are 
primarily quasars. In Figure [3] we display a few example 
BOSS spectra of our SN host galaxies, spanning the full 
redshift range of SDSS-II SNe. Though there is a wide 
range of S/N in these spectra, spectral features (partic- 
ularly the 4000 A break and emission lines), which allow 
us to measure the galaxy redshift, are clearly visible. 

3. CREATING A PHOTOMETRIC HUBBLE DIAGRAM 

In this section, we describe our construction of a 
photometrically-classified SN la Hubble diagram. We 
provide details of the SN simulations used to determine 



the optimal selection criteria (Section 13. ip and apply 
these to the data (Section |3~2|). We describe the resulting 
Hubble diagram in Section [3T3l 

3.1. Selection Criteria 

Here we define, using simulations, the criteria by 
which we construct our photometric SN la Hubble dia- 
gram. Our primary focus is to minimize the contamina- 
tion in our sample from non-la SNe; this is one of the 
major concerns associated with measuring cosmological 
parameters with photometrically-classified SNe, and has 
the potential to introduce significant systematic errors 
in the cosmological analysis (Section |6]) . Given the large 
number of SN candidates included in our data-set, the 
conservative classification criteria that we seek will still 
result in statistical errors on our cosmology smaller than 
those due to systematics (Appendix |Pj) . By focusing pri- 
marily on the purity of our sample we have sacrificed its 
overall completeness, and therefore we caution the reader 
against using this particular sample for analyses that re- 
quire high completeness (e.g., SN rate measurements). 

There are three types of criteria that we apply to 
our sample. First we apply "light-curve quality" cuts 
which are criteria that ensure the SN light-curves (real 
or simulated) are of sufficient quality to be well-fit with 
SALT2 (see Appendix IB"! for details) to give accurate dis- 
tance moduli and provide meaningful typing constraints. 
Next we consider the optimal values for the PSNID clas- 
sification parameters (see Appendix [A] for details), which 
will differ from those used in Sll since we are both using 
an improved version of PSNID and wish to increase the 
recovered sample purity. Finally, we can further improve 
the purity of our sample by applying "color and stretch" 
criteria based on the derived SALT2 parameters for each 
light-curve and our knowledge of the likely acceptable 
range for these properties for SNe la. We outline each of 
these types of criteria below, following a brief overview 
of the simulations and Figure-of-Merit (FoM) we use to 
guide our criteria selection. In Table [IJ we show the ef- 
fect of each selection criterion as applied in the following 
sections. 

3.1.1. Supernova Simulations 

To test the purity and completeness of our sample, 
and help define the best selection criteria, we use the 
publicly available simulations of the SDSS-II SN Survey 
created as part of the "Photometric SN Classifier Chal- 
lenge" of KlOa. These simulations were made available 
via the SN challenge websit and were p roduced using 
the SNANA software (jKessler et al.ll2009al) . as described 
in K 10b. We decided to use these simulations to test our 
completeness and purity as they have been well-tested by 
many researchers, and were designed specifically for test- 
ing photometric classification of SNc. They also provide 
an accurate description of the conditions under which 
our data were acquired. 

In detail, the simulations have ten times the number 
of SNe as the full three-year SDSS-II SN Survey. They 
are based on realistic weather, seeing, and photometric 
zero-point variations, and have a realistic mixture of dif- 
ferent SN types out to a redshift limit of z = 0.45. The 

2 http://sdssdp62.fnal.gov/sdsssn/SIMGEN.PUBLIC 
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Fig. 3. — Galaxy spectra from BOSS. The top left panel is a low redshift galaxy with a high S/N continuum, while the top right panel 
shows a galaxy with virtually no continuum, but several clear emission lines. Both bottom panels are at the high end of our redshift range 
(z>0.4). The black lines show the data, the green is the best-fit eigenspectrum spectra and the errors are in red (except masked points, 
which are set to zero). 



number of SNe la created in these simulations is based on 
the observed SN la redshift depende nce from the SDSS-I I 
SN Survey first-year rates analysis (jDildav et al.l l2008h . 
and the non-la SNe redshift dependence is based on 
the c ore-collapse rate analysis from SNLS ([Bazin et al.l 
f2009h . The core-collapse contribution has been inten- 
tionally overestimated in the simulations in order to in- 
crease the statistics of non-la SNe that are misidentificd 
as SN la, rather than underestimate the amount of non- 
la contamination. 

For the analysis presented herein we only use the 
SALT2 simulated SN la light-curves (ugriz) in this pub- 
lic data-set (the simulations also included MLCS gener- 
ated light-curves). Since this is only half of the gener- 
ated SNe la in the simulation, we only include half of 
the non-la SNe to keep the ratio of SNe la to non-la SNe 
correct. The SNe la light-curves use the SALT2 stan- 
dardization parameters of a = 0.11 and /3 = 3.2, with 
an assumed intrinsic dispersion of ai n t = 0.12 mag. The 
intrinsic dispersion is included in the simulations of the 
SNe la by adding random color variations for each SN in 
each passband, drawn from a Gaussian distribution with 
a m = 0.09 mag, applied coherently to all SN epochs. 
The non-la SNe in the simulations are based on 41 well- 
measured spectroscopically-confirmed non-la SNe tem- 
plates. A flat ACDM cosmology with fl m = 0.3 and 
Ho = 70 km s _1 Mpc -1 was also assumed. 

We note that the simulations accurately model the 
SDSS-II SN Survey software detection pipeline, and thus 
any SNe that were too faint to be detected by the SDSS 
survey would not be included in the simulations. The 
simulations can also model spectroscopic selection, which 
we do not make use of in our main analysis but is included 
where we make checks against our spectroscopically- 



confirmed subsample. 

3.1.2. Defining the Figure-of-Merit 

As we know the true types of SNe in the simula- 
tions, we can estimate the efficiency and purity of our 
photometric classifications as a function of our selection 
criteria. We use the definition of photometric typing ef- 
ficiency ei a from Sll, 



Aftrue 



(i) 



where A^ UT is the total number of true SNe la in the 
simulations that pass the light-curve quality cuts (see 
Section [5X3)1, and N^ uc is the subsample of iVg UT that 
are classified correctly (itself a function of additional 
cuts; see Tabic [TJ. As noted in Sll, this equation mea- 
sures the efficiency of classification only for well-observed 
SNe la, not the total efficiency of identifying all SNe la 
in our simulated data. From here on in we refer to the 
parameter defined in Equation [1] as simply our efficiency. 
We additionally define the weighted purity as in Sll, 



'/In 



Wj false 



^yfalsc 



(2) 



where N(^ lsc is the number of non-la SNe incorrectly clas- 
sified as SNe la and W{^ lsc weights the contribution of 
misclassifications to the overall purity. This definition 
has the usual meaning of purity for Wj* lsc = 1 , and for a 
given amount of contamination by non-la SNe a higher 
(lower) value of results in a lower (higher) value 

of TJi a . 

The purity and efficiency of the photometric classi- 
fication can be combined to form a FoM. As defined in 
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Sll and KlOb, the FoM is simply the product of Eqns[T] 
andU 

7W"true j\rtruc 
p„J/_ la iv Ia /o\ 

jV CUT (N tmc + w {alsc N ialsc ) ' 
la \ la la la / 

The FoM in Eqn [3] does not encapsulate information on 
the cosmology constraints, but rather is a simple metric 
that describes the broad merits of a classifier. Never- 
theless, we aim to create a classification that optimizes 
this FoM with a suitably chosen weig hting factor (Wff se ) 
that represents our previously stated choice of prioritiz- 
ing purity over efficiency. The ideal choice for W{^ sc is a 
complicated function of the contaminating objects, sen- 
sitive to the redshift and magnitude distribution of each 
subtype. This issue is not investigated in detail here; in- 
stead, we tested the effects of several different weighting 
values greater than one and empirically determined that 
Wja = 5 is the best choice, as it produced the sharpest 
peak in the FoM in Figures HI [5] and [7] We therefore 
use this value in all subsequent analyses. For further 
discussion of the importance of purity in SN la sam ples 
for cosmological an alyses, see lBernstein et ahl (|2011l ) and 
iGiergo et al.l (|2012h . who discuss this issue in the context 
of the Dark Energy Survey. We note that the efficiency, 
purity and FoM are only calculated after we have ap- 
plied our light-curve quality cuts, as this is our baseline 
for defining SNc la that arc potentially useful for cosmo- 
logical constraints. 

3.1.3. Light-curve Quality Cuts 

We begin our analysis of the simulations by apply- 
ing data quality cuts to all light-curves, removing SNc 
that have insufficient epochs to provide any useful mea- 
surement. Defining t as the rest-frame epoch (in days) 
of each SN relative to peak (determined from the best-fit 
PSNID SN la model), we require at least one epoch of 
photometry near peak at — 5 < t < +5 and at least one 
additional epoch at t > 15, as in Sll. However, we do 
not apply the S/N criteria outlined in Sll (S/N > 5 in 
at least two of the gri bands), as this could remove many 
high-redshift or under-luminous SNe that may be iden- 
tified through photometric classification. We note that 
an implicit S/N limit does in fact exist, as the difference 
imaging software (sdssdiff ) in the SDSS-II SN Survey 
requires multiple detections at S/N > 3 for a n object to 
be la belled a SN candidate in the first place (jSako et al.1 
[2001 . Finally, we check again that none of our SN can- 
didates were detected in more than one of our SDSS-II 
SN search season. 

Since we have no hard S/N limit, it is not strictly 
true that our light-curve quality cuts remove all objects 
that are incapable of providing useful cosmological con- 
straints from their light curves. The simplicity of our cri- 
teria is an attempt to balance the necessity of sufficiently 
useful data with the desire to be unbiased against faint 
objects. 

3.1.4. PSNID Criteria 

We run PSNID (described in Appendix [5} on all 
simulated light-curves that pass our light-curve quality 
cuts (Section 3.1.3), placing flat priors on Ay (the host- 
galaxy extinction), T max (the time of peak brightness), 
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Fig. 4. — The efficiency (blue dashed line), purity (red dot-dashed 
line), and FoM (green solid line) for the simulated sample as a 
function of the position of the PSNID Pi a probability cut. We plot 
the true purity (W^ lsc = 1), and only change the weighting factor 
to W^ lsc = 5 in the FoM. 

and A?Tii5 (the stretch parameter). We use the true red- 
shift of each SN as a prior, with an uncertainty on z of 
the measured error. 

We investigate PSNID criteria for removing non-la 
SNe using our simulations, examining their effect on the 
efficiency, purity, and FoM. First we optimize the cut on 
the PSNID probability of being a SN la (P Ta ). In FigureH 
we show the efficiency, purity, and FoM as a function of 
Pi a . Due to the general behavior of PSNID, which tends 
to cluster values of Pi a around zero or one (demonstrated 
in Figure 7 of Sll), these functions are relatively flat. As 
Pia does not provide much discriminating power beyond 
these extreme values, the FoM has little sensitivity to the 
Pi a cut value. We thus require Pi a > Pibc and Pi a > P\\ 
for inclusion in our SN la sample, which combines high 
efficiency with modest purity. For reference, this is a less 
constrained criterion than Sll, which adopted P[ a >0.9. 

We also use the reduced chi-squared (Xr) °f the 
best-fit PSNID model as a discriminator of non-la SNc, 
and determine the optimal value for this cut after the 
PSNID cut (P a > Pibc and P[ a > P n ) is applied. In 
Figurc[5]one can see a clear peak in the FoM at Xr — 1-2, 
close to where the purity and efficiency curves cross. 
Thus wc define a PSNID classified SNe la to be an object 
with P a > Pibc, Pa > Pi, and xl < 1-2 . 

The value of our Xr cut differs significantly from that 
in Sll, which can be seen in their Figure 10 to be located 
at a broad maximum of Xr — 1-8- As noted in Appendix 
A, the version of PSNID used herein differs from that of 
Sll (which had larger model uncertainties), resulting in 
the optimal Xr cut being smaller in this work. 

3.1.5. SALT2 Criteria 

Despite our optimization of the PSNID classification 
criteria, over 25% of our photometric SN la sample re- 
mains non-la SNe (see Figure [5] and Table [1]). The purity 
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Fig. 5. — The efficiency (blue dashed line), purity (red dot-dashed 
line), and FoM (green solid) for the simulated sample as a function 
of Xr> after the PSNID cut (Pi a > P Ibc and Pi a > Pa) was ap- 
plied . We plot the true purity (W^ lae = 1), and only change the 
weighting factor to \yt^ se = 5 in the FoM. Our cut at Xr = 1-2 is 
denoted by the verticaf black line. 



Fig. 6. — The X\ and color (c) distributions for the simulated 
SNc remaining in our photometric sample after PSNID cuts with 
SNe la as black points and non-la SNe as blue cross symbols. The 
purple ellipse defines the area we keep for the photometrically- 
classified sample of SNe la. The cutoff at —5 and 5 in X\ are hard 
limits set by SALT2. 



of this sample is unsuitable for the cosmological analysis 
discussed in Section |H which is the goal of this paper. 
We thus run the SALTS light -curve fitter (described in 
detail in Appendix [B]) on all of our SNe (simulated and 
data) to obtain the best-fit light-curve parameters, and 
explore the effects of applying additional selection cri- 
teria to these parameters to further differentiate SNe la 
from non-la SNe. 

In Figure |6]wc show the distribution of the measured 
SALT2 parameters color (c) and "shape" (Xi) for all SNc 
remaining in our photometric SN la sample. By defini- 
tion, SNe la form a well-defined cluster of points centered 
on zero in this par ameter space, while no n-la SNe are 
more scattered. In iLampeitl et al.l (|2010al) independent 
limits were placed on X\ and c, but it is c lear that an 
ellips oidal cut (similar to the circular cut of iBazin et al.l 
1201 1| ) would yield a higher SN la purity for a given ef- 
ficiency. We determine the optimal lengths of the semi- 
major (clx-l) an d semi-minor (a c ) axes of this ellipse in 
Figure [3 which shows the efficiency, purity, and FoM as 
a function of the ellipsoidal parameters. For an ellipse 
centered at (xi,c) = (0,0), the FoM shows a clear peak 
at a Xl = 3 and a c = 0.25. We use this ellipse, shown 
in Figure [SJ to further distinguish SNe la from non-la 
SNe. As can be seen in Table [TJ this procedure removes 
~ 70% of our contaminating SNe at the expense of re- 
jecting ~ 20% of the SNe la. We investigated allowing 
the centre of the ellipse to vary, but found this did not 
significantly improve the purity or efficiency. 

3.1.6. Color-magnitude Cut 

3 We use only SALT2 throughout this paper. As we are inter- 
ested in the utility of photometric samples of SNe la for cosmol- 
ogy, rather than presenting definitive cosmological constraints, we 
forego a detailed comparison of light-curve fitting algorithms. 



IBazin et al.l (|201lD recently showed that a cut in 
observed color-magnitude space significantly improves 
SN la sample purity by removing core-collapse SNe con- 
taminants. We investigated the effects of such a color- 
magnitude cut on our simulations using the gri peak 
magnitudes from the best-fit SALT2 model for each SN. 
We found the most effective color-magnitude cut is in the 
g — r versus g— band model magnitude plane. The opti- 
mal orientation of this color-magnitude cut is described 
by g — r < 0.3 x (g — 21.2), and is shown in Figure [H 
This diagnostic cleanly removes a population of non-la 
SNe in the simulations which are too "blue" at the given 
magnitude to be SNe la. We found that the application 
of additional color-magnitude constraints, e.g., using dif- 
ferent filters, did not produce a significant reduction of 
our contamination beyond the first color-magnitude cut. 

3.1.7. Overall Contamination in the Simulations 

In Table [T] we break down the effects of our selection 
criteria on our SN sample. The first row of Table [1] lists 
the number of SNe of each type that are included in the 
simulation; successive rows contain the number of SNe 
that remain in our SN la sample after each of the criteria 
(as detailed in the previous sections) arc applied. We 
provide the contamination (1— purity, assuming W^ lse = 
1), efficiency, and FoM (assuming W^ lsc = 5) after the 
application of each selection cut. These quantities are 
defined using the number of SNe la that pass our light- 
curve quality cuts as the total number of SNe la, i.e., 
Nf a VT in Eqn 1. 

Application of both the PSNID and SALT2 criteria 
significantly improves the purity of our sample, result- 
ing in a SN la sample with purity > 90% and efficiency 
> 70%. The FoM has also increased significantly with 
these cuts, reflecting the weighting we have applied to pu- 
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Fig. 7. — The efficiency (top), purity (middle), and FoM (bottom) 
plots from SN simulations, for changing the semi-major axis (oxj 
and semi-minor axis (a c ) of the ellipse in Figure[6] 

rity. Interestingly, the inclusion of the color-magnitude 
constraints (Figure [5]) significantly improves our purity 
with almost no effect on our efficiency. 

Figure [5] is the Hubble diagram for our simulated 
photometrically-classified SN la sample. The plotted er- 
rors are a combination of the uncertainties on the SALT2 
light-curve parameters, the redshift and an assumed in- 
trinsic dispersion of Oi n t = 0.12 magnitudes. The SALT2 




17 18 19 20 21 22 23 24 



r 

Fig. 8. — The SN g — r color as a function of peak SN g— band 
magnitude derived from the best-fit SALT2 model to our simula- 
tions. The black dots are SNe la and the blue cross symbols are 
non-la SNe. The purple line is our color-magnitude cut that best 
separates non-la SNe contaminants from our sample. 



SN la parameters are fixed to the input simulation values 
of a = 0.11 and /3 = 3.2. 

We show in blue in Figure [9] the 106 misclassificd 
SNe that have passed all of our selection cuts. This final 
simulated photometric SN la sample has a contamina- 
tion of 3.9% and an efficiency of 70.8%. Our purity is 
higher than that for any of the participating methods in 
the "Photometric SN Classifier Challenge" (KlOb), al- 
though it is necessary to note that we have explicitly 
placed a higher priority on purity than was the stated 
goal in KlOb, and their analyses used a DES-like set of 
simulations that extend to higher rcdshifts than in our 
sample. As expected, the source of our contamination 
is primarily SNe Ib/c, as their light-curves most closely 
resemble those of a SNe la. 

Finally, in Figure [TQl we show the expected purity 
and efficiency of our sample (based on simulations) as a 
function of redshift. Within the errors, we see no signifi- 
cant redshift dependence in our purity, but we do witness 
a significant fall in our efficiency at higher redshifts. This 
is expected, as we have prioritised purity over efficiency 
for our cosmological study. Larger simulations (see Sec- 
tion 7.2) are required to probe more subtle effects with 
redshift. 



3.2. Application of selection method to the data 

We now apply our analysis and the selection criteria 
defined using our SN simulations to our SDSS-II BOSS 
SN sample. In Table [2J we provide a break-down of the 
number of SNe classified at each stage of our selection (as 
discussed in Section l3~Tj) . The two PSNID-based criteria 
remove sa 45% of the sample that remain after the data 
quality cuts we applied (comparable to the 48% seen in 
our simulations), resulting in 1443 objects in our SN la 
sample at this stage of analysis. 
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Cut 


Classified SNc la 


Contam- 


SN la 


FoM 




Total 


SNe la 


non-la SNe 


ination 


Efficiency 


(W=5) 


Number in simulation 


12203 


5018 


7185 








Light-curve quality 


9186 


3734 


5452 


59.4% 


100% 


12.05 


Pla > Pnon-Ia (PSNID) 


6354 


3701 


2653 


41.8% 


99.1% 


21.6 


xl <l-2 (PSNID) 


4737 


3420 


1317 


27.8% 


91.6% 


31.3 


X\ and c cut (SALT2) 


2918 


2675 


243 


8.3% 


71.6% 


49.2 


Color— magnitude criteria (SALT2) 


2750 


2644 


106 


3.9% 


70.8% 


59.0 



TABLE 1 

Number of SNe in our simulated sample as a function of the selection criteria applied to the data. In each row we show 
the cumulative effect of all the previous criteria on the contamination, efficiency and fom (assuming w^ laa = 5). 



Selection Criteria 


Removed SNe 


Kept SNe 


Spec SNe la 


Spec non-la SNe 


Accurate BOSS Redshifts 




3500 


329 


59 


Light— curve quality 


874 


2626 


249 


24 


Pla > Pnon-Ia (PSNID) 


579 


2047 


247 


6 


xl <l-2 (PSNID) 


604 


1443 


239 


2 


X\ and c cut (SALT2) 


634 


809 


209 





Color— magnitude criteria (SALT2) 


54 


755 


209 





Correct host galaxy 


3 


752 


208 






TABLE 2 

The cumulative effect of applying each selection criteria to OUR DATA, LEADING TO A FINAL SAMPLE OF 
PHOTOMETRICALLY-CLASSIFIED SNE IA (752). THE RIGHT-HAND COLUMN SHOWS THE EFFECT OF THESE CRITERIA ON KNOWN 

(SPECTROSCOPICALLY-CONFIRMED) SNe IA IN OUR SAMPLE. 




Redshift Redshift 



Fig. 9. — Hubble diagram of our simulated SNe la that pass all of 
our selection criteria. This includes 106 non-la SNc that have been 
misclassificd (blue symbols) and 2644 correctly classified SNe la 
(black symbols). The redshift limit of z = 0.45 is artificial and set 
by the original limit in the public SNANA simulations. The plotted 
errors are a combination of the uncertainties on the SALT2 light- 
curve parameters and the redshift. The SALT2 SN la parameters 
have been fixed to the input simulation values of a = 0.11, /3 = 
3.2, and M=29.8. The bottom panel shows the Hubble residuals 
assuming the input cosmology. 

We use SALT2 (Appendix E}, applied to the SDSS 
SMP griz light-curve data, to measure the light-curve 
parameters and distance moduli for all 3500 (2382 + 941 
+ 177) SN candidates with host galaxy redshifts (see 
Section 12. 3[) . Figure [11] shows the distribution of our 
PSNID-classificd SN la sample in Xi-color space and the 
ellipsoidal criteria we derived from simulations. This cut 



Fig. 10. — The simulated efficiency (blue) and purity (red) for 
our final photometrically-classified sample. The error bars are de- 
termined from the propagation of errors using Eqns [l] and [2] 

removes 634 SN candidates (44%) from our photomet- 
ric sample. This is a higher percentage than predicted 
by simulations in Section I3TT1 (39%), which could be at- 
tributed to the fact that non-SN contaminating sources 
(e.g., quasars, which are known to be in our BOSS tar- 
get list) are not modeled in our simulations. Many of the 
discarded SN candidates have poorly-fit templates, with 
Xi = ±5, where SALT2 is driven to the extremes of its 
self-imposed (i.e., hard-coded) X\ range. This constraint 
also removes many "red" objects with large c values, 
which may include highly reddened SNc la. However, we 
would prefer to exclude these red SNe la from our cos- 
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Fig. 11.— The Xi and c distributions for PSNID-classificd 
SNe la from our SDSS-II SN candidate samples. The blue (red) 
cross symbols denote the sub-set of these candidates that have been 
spectroscopically confirmed as SNe la (non-la SNe). The purple el- 
lipse is our SNe la boundary taken from Figure [6) 
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Fig. 12. — The g — r color as a function of g-band magnitude for 
our SDSS-II data. The blue cross symbols denote the sub-set of 
these candidates that have spectroscopic confirmation as SNe la. 
The location of the purple line is our SN la boundary taken from 
Figure [8] 



mology analysis, as lFolev et al.l (|2010bf ) showed they po- 
tentially bias cosmological constraints. Our SALT2 cri- 
teria leave us with 809 photometrically-classified SNe la 
candidates. 

In Figure [12] we show the application of the color- 
magnitude cut in the g — r versus g— band model mag- 
nitude plane, which removes a sample of "blue" SN can- 
didates. Both in the data and the simulations this cut 
removes ~ 7% of the sample that passed the SALT2 cut, 
leaving a photometrically-classified sample of 755 SNe la 
candidates. The agreement between the simulations and 
data is reassuring, and it should be noted that these sim- 
ulations were made prior to the construction of our BOSS 
SDSS SN sample, so there has been no fine-tuning of the 
simulations to match our BOSS sample. 

A potential source of error that is not included in 
our simulations is the misidentification of the host galaxy 
of our SN candidates, as this can result in an incorrect 
classification due to an incorrect rcdshift prior. We have 
identified three cases in our sample where we have evi- 
dence, described in Appendix [C] that either the galaxy 
observed by our BOSS program is not associated with 
the SN to which it was assigned, or the derived BOSS 
rcdshift is erroneous. We remove these three SN candi- 
dates from our sample. 

Our final sample contains 752 photometrically- 
classified SNe la, of which 208 have been spectroscopi- 
cally confirmed as SNe la as part of the original SDSS-II 
SN Survey (Tabled]). 

3.3. Tests of our photometric sample 

We present here a basic examination of our final 
sample of 752 photometrically-classified SNe la before 
embarking on a cosmological analysis (Section [5]). First, 
we study the effect of our selection criteria on the sub- 
set of existing spectroscopically-confirmed SNe la in our 



sample. Unlike Sll, we have not used this subset of 
known SNe to refine our selection criteria as there are 
concerns about potential bias in this sample (see Sec- 
tion 12. ip . As can be seen in Table [5] we start with 329 
spectroscopically-confirmed SNe la that have had host- 
galaxy redshifts obtained by BOSS or SDSS. Of these, 
249 SNe la passed our data-quality criteria applied in 
Section 13.1.31 After applying all of our additional se- 
lection criteria we are left with 208 spectroscopically- 
confirmed SNe la, resulting in an efficiency of 84% (208 
/249) in classifying our SDSS-II spectroscopic SN la sam- 
ple. 

This efficiency is higher than the predicted value 
from our simulations (71.6%; see Table [JJ. This differ- 
ence can be mainly attributed to the spectroscopic sam- 
ple probing a lower redshift range than the photometric 
sample. In Figure [13] we show the redshift distribution 
for our full sample of photometrically-classified SNe la 
(black) compared to the subsample of spectroscopically- 
confirmed SNe la (blue). While the spectroscopic sam- 
ple peaks at z ~ 0.2 and drops to zero by z > 0.4, the 
photometrically-classified SNe la extend out to z ~ 0.55, 
with a median redshift of z = 0.30. This explanation is 
checked by studying the "spectroscopically-confirmed" 
subset of SNe la provided as part of the public SNANA 
simulations discussed in Section 3.1.1. The efficiency of 
photometrically classifying this simulated spectroscopic 
SNe la subset is 83.2%, in reassuring agreement with the 
84% efficiency we find for the spectroscopic SNe la subset 
in our data. 

We reiterate that we do not include in our final 
sample any of the spectroscopically confirmed SNe la 
that were removed during our photometric classification 
procedure. Our intention is to construct a sample of 
SNe la based purely on their photometric properties and 
host galaxy redshifts, and examine the cosmological con- 
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Fig. 13. — The redshift distribution of our 752 photometrically- 
classified SNe la (black), and the subset of 208 that have spectro- 
scopic classification (blue). 

straining power of such a sample. This procedure mimics 
the challenges facing the next generation of SN surveys. 

All spectroscopically-confirmed non-la SNe in the 
original SDSS-II SN sample were removed by our se- 
lection criteria, resulting in 100% purity for the spec- 
troscopic subset of our photometric sample. Of the 59 
known non-la SNe with host-galaxy redshifts obtained by 
BOSS, only two were misclassified as SNe la by PSNID, 
and both were subsequently removed by the X\ ~ c cri- 
teria. There are two factors that lend themselves to this 
subset achieving higher purity than would otherwise be 
expected from our simulations. First, the SDSS-II spec- 
troscopic program was intentionally biased against tar- 
geting non-la SNe; in our BOSS sample there is only 
one spectroscopically-confirmed non-la SN for every six 
SNe la, whereas nearly 60% of the SNe in our simulations 
are non-la. Second, PSNID was tested, and optimized, 
on the SDSS-II spectroscopic SNe, so its efficiency at 
identifying non-la SNe from this sample is unsurprising. 

In Figurc[14l we show the SALT Xi values as a func- 
tion of redshift for the photometrically classified SNe la 
sample. We also show the mean Xi in redshift bins and 
the best-fit linear relation to the binned data. We see no 
change in the distribution of X\ values with redshift, as 
the best-fit slope is consistent with zero within the errors 
on the fit. In Figure [15J we show the SALT c (color) val- 
ues as a function of redshift, again with their mean val- 
ues, binned by redshift, and the best-fit linear relation. 
A correlation is evident in this case, with higher redshift 
SNe skewed towards bluer (negative) colors. However, 
this trend is driven by SNe at z > 0.4, and if we limit 
the sample to redshifts below this value we find the best- 
fit slope is consistent with zero within errors. We may 
be seeing evidence of a color-dependent Malquist bias, 
i.e., bluer SNe la are brighter and thus easier to detect 
at higher redshifts (more so than SNe with higher X\ 
values) . 

In Figure[l6]we present two Hubble diagrams for our 
photometrically-classified sample of SNe la. On the left 
is the sample before our color-magnitude cut is applied, 
while on the right we show our sample after making this 
cut. This criterion removes 54 SNe la candidates, which 
are clearly offset (too faint) from the main Hubble se- 
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Fig. 14. — The SALT2 X\ parameter as a function of redshift for 
the photometrically-classified SNe la (black dots) and the subset 
of spectroscopically-confirmed SNe la (blue crosses). The green 
points show the mean X\ (and the error on the mean) in bins of 
redshift. The red solid line is the best-fit linear relation to the 
average X\, and the red dashed lines are the error on the fit. 

quence and thus likely to be non-la SNe (as seen in our 
simulations). These removed candidates are also clus- 
tered around z w 0.2, again consistent with our simula- 
tions. 

There are still ten possible outliers, around z « 0.2, 
to the main Hubble diagram in Figure [TO] (right-hand 
side). We have studied the photometric data for these 
HR outliers individually (including visually inspecting 
the images) and can find no obvious reasons to remove 
them from our sample. These outliers could be Type 
Ibc SNe, which possess similar colors and light-curve 
shapes in the SDSS filters at these redshifts, or possibly 
odd SN e la like PTFlOops (|Maguire et al.ll20l"TI) and SN 
2006bt (jFolev et al.ll2010al) . However, this small amount 
of contamination does not bias our cosmological results, 
as shown in Section I6.1[ where we apply at 3a clipping 
to the Hubble diagram and find consistent cosmological 
results. 

4. SOURCES OF BIAS 

Before we undergo a detailed cosmological analysis 
of our sample (Section|H|), we investigate possible sources 
of bias in our photometrically-classified sample. 

In Figure [T7] we show the residuals about the Hub- 
ble diagram from our simulations (Section [XT}; the Hub- 
ble residual (HR) is the difference between the observed 
and input distance modulus for each SN la. We have 
subtracted the true distance modulus, based on the in- 
put cosmology, from the measured distance modulus for 
each SN in the simulation, and binned the results in red- 
shift. The left panels of Figure [T7] show the SNe la 
photometrically-classified using the methodology out- 
lined in Section (3[ while the right-hand panels shows 
residuals for the true SNe la in the simulations. Each 
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Fig. 15. — The SALT2 color parameter as a function of redshift 
for the photometrically-classified SNc la (black dots) and the sub- 
set of spectroscopically confirmcd-SNe la (blue crosses). The green 
points show the mean color (and the error on the mean) in bins 
of redshift. The red solid line is the best-fit linear relation to the 
average color, and the red dashed lines are the error on the fit. 



row demonstrates the effect of adding one of our selec- 
tion criteria, described in Section |3~T1 In the absence of 
any systematic bias, we would expect these residuals to 
be scattered about zero (blue line). 

There are two sources of bias whose effects are eas- 
ily seen in Figure [17] First, after making only the data- 
quality and PSNID cuts, our SN sample (second row, left 
panel) shows a strong bias towards positive HRs (under- 
luminous objects) at z < 0.2. We have shown in Tabled] 
that this sample has high contamination, i.e., at z < 0.25 
non-la SNe, though fainter than SNe la, are still bright 
enough to be detected as SN candidates in the SDSS 
SN Survey. At higher redshifts core collapse SNe are 
too faint to be observed, and thus their contamination is 
naturally curtailed. This bias at low-z, due to contam- 
ination, is effectively eliminated once the Xi-color cut 
is introduced. After the final color-magnitude cut (bot- 
tom panel) the HRs are consistent with zero bias at low 
redshifts. 

At high redshifts (z > 0.25) the HRs in all pan- 
els of Figure [T7] are biased low; i.e., SNe are, on av- 
erage, brighter than expected given the input cosmol- 
ogy. This selection effect is a comb i nation of the clas- 
sic M almquist bias (|Malmquistlll936l : IWood-Vasev et al.1 
I2007h and other SALT2-related selection effects. We 
have only analysed objects that would be classified as 
SN candidates in the SDSS-II SN Survey (survey pipeline 
cuts and post-survey cuts are applied within the simu- 
lations), so under- luminous objects at high-z are never 
output from the simulations. Additionally, though we 
place no implicit S/N cut on our data, low S/N SNe 
that are poorly fit by SALT2 are more likely to fail the 
X\ — c cut, and our original data-quality cut, and are 
thus preferentially excluded (see Section [STTj) . 

For clarity, we show in Figure [18] the differ- 



ence between the recovered distance moduli for the 
photometrically-classified SNe la and the true SNe la 
as a function of redshift, i.e., the bottom left-hand panel 
in Figure [17] minus the bottom right panel in Figure 1171 
This removes the Malmquist bias effect, allowing us to 
see the residual contamination bias as a function of red- 
shift. As one would infer from Figure [TU] the redshift- 
dependent bias of our HRs due to contamination is seen 
in Figure [18] to be a small eff ect. We find, using Akaike 
information criterion (AIC; lAkaikel Il97l to compare 
the two models, that a linear fit is favored over a redshift 
independent model, with a slope of -0.021 best describ- 
ing the linear bias in the measured redshift range. This 
potential redshift dependence deviates from null by only 
—0.01 mag out to z = 0.5. Thus the effect of the contam- 
ination is subdominant to the Malmquist bias, which has 
modeled uncertainties that are larger than the possible 
redshift dependence of the contamination bias. 

The combination of these selection effects (the clas- 
sical Malmquist bias and SALT2 effects) are by far the 
biggest source of bias affecting photometric SN sample 
if left uncorrected, which is discussed in detail in Sec- 
tion |4J] However, we note that such effects, especially 
the classical Malmquist bias, are also present in spectro- 
scopic samples where they are more difficult to correct 
for, as this bias will depend on the details of the spectr o- 
scopic program (|Sako et al.ll200llKessler et al.ll2009aD . 

4.1. Selection effects 

In this section we correct for the known selection 
bias at z > 0.25 seen in FigurefTTl This bias is a combina- 
tion of the classical Malmquist bias and SALT2-related 
effects. For a magnitude limited SNe survey, the classical 
Malmquist bias results for a given X\ and color in the 
preferential detection of SNe that appear brighter due 
fluctuations caused by Poisson noise and intrinsic scat- 
ter. We do not attempt to study these different effects 
separately in this paper and only provide a correction for 
their combined effect on the distance moduli of our pho- 
tometric SNe la sample. For this reason, we will from 
hereon refer to this combined selection effect as just the 
"Malmquist bias" , but ask the reader to recognise that it 
is a combination of magnitude (or S/N) effects, of which 
the classical Malmquist bias is likely the most important. 

Before we present the details of our correction to the 
Malmquist bias, we first demonstrate the importance of 
this effect in our cosmological analysis. We compute the 
equation-of-state of dark energy (w) for the entire sam- 
ple of true SNe la that pass data-quality cuts (top right 
panel of Figure [TT)) using the publicly ava ilable software 
package CosmoMC (]Lewis fc Bridldl2002T) . including pri- 
ors on the cosmological parameters (we detail our further 
usage of CosmoMC in Section [6J. We find the best- 
fit cosmology for this sample to be w — —0.90 ± 0.05, 
which is inconsistent with the input cosmology of w = — 1 
(though only at the two-sigma level), and may indicate 
a bias. We stress that this sample is, by definition, 100% 
pure and 100% efficient, yet still produces biased cosmo- 
logical parameters. 

Therefore, to investigate the selection effects bias, 
we have undertaken a new set of SNANA simulations 
that span a wider range in redshift than those used in 
our analysis of the sample contamination (Scction l3.1.1[) . 
We also increase the number of SNc in the simulation, as 
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Fig. 16. — Hubble diagram for our photometrically-classified SNe la sample. Blue points are the subsample of spectroscopically-confirmcd 
SNe la, while black points are new SNe la. These Hubble diagrams are created using the best-fit a and /3 from the full cosmological fit 
(Section |B), and include the intrinsic dispersion in the uncertainties. Our sample before (left) and after (right) the color— magnitude criteria 
are applied demonstrates the utility of this criterion. The SALT2 SN la parameters have been set to the best-fit values of a and /3 from the 
cosmological fit, and as our cosmological fit analytically marginalizes over M we use the same value here as in the simulations (M=29.8). 
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Fig. 17. — The weighted mean and uncertainty for Hubble resid- 
uals in our simulations in bins of Az = 0.02. The left-hand panels 
show the photometrically-classified SNe la, while the right-hand 
panels are for true SNe la. Each row shows the cumulative effec t 
of applying the various selection criteria discussed in Section 13.11 



Fig. 18. — The difference between the weighted mean of the Hub- 
ble residuals from the photometrically classified SNe la and the 
true SNe la, i.e., the difference between the bottom two panels in 
Figure 17. The error bars account for the correlations between the 
photometrically classified SNe la and the true SNe la. 



the Malmquist bias is small compared to the size of the 



error-bars and requires sufficient S/N to characterize it. 



16 



0.05 




0.2 0.3 
Redshift 



Fig. 19. — The difference between the observed (/i b s ) an( i ex- 
pected (/icosmo) distance modut us as a function of redshift for the 
simulations described in Section 14.11 We show the best-fit expo- 
nential function to this data (blue line) and the one-sigma errors 
on the fit (red dashed line). 



In det ail, we use SNANA (version v9_97; iKessler et al.1 
l2009a| ) to create thirty thousand SNe la over the redshift 
range 0.01 < z < 0.5. For consistency, the SALT2 SN 
parameters and the assumed cosmology are the same as 
those used for the public SNANA simulations described 
in Section 13.1.11 As we are characterizing an effect for 
SNe la, we do not include core-collapse SNe in the sim- 
ulation. 

In Figure \T9\ we show the weighted mean differ- 
ence between the "observed" (nobs) and true (/i c05m o) 
distance modulus for our simulated SNe la after apply- 
ing our photometric-classification criteria to the datey. 
We use Eqn IBll to calculate Mobs, fixing a and j3 to the 
values used in the simulatioro There is little bias out 
to z ~ 0.3, but at higher redshifts a significant offset 
appears, growing with increasing redshift and becoming 
> 0.1 magnitude over- luminous. 

We also show in Figure [19] the best-fit analytical 
function to the Malmquist bias, which is an exponential 
of the form, 



/^corr(^) 



ae bz + c, 



(4) 



where a = -0.004 ± 0.001, b = 7.26 ± 0.31, and c = 
0.004 ±0.006. We include this Malmquist bias correction 
in our distance modulus calculation (Eqn lBTj) . resulting 

4 We apply all our data cuts to the simulated data except for the 
PSNID cuts, which require significant computational time for such 
a large suite of simulations. We have confirmed that excluding this 
cut does not effect our results using the public SNANA simulations 
since we are only interested in using true SNe la for the Malmquist 
bias correction. Only an additional 7% of SNe la in these simu- 
lations are removed by this cut. We have also demonstrated that 
the shape of the Malmquist bias is no different whether or not the 
PSNID cut is included. 

5 Allowing a and /3 to vary during the computation of fi does 
not significantly change our Malmquist bias correction, but does 
increase the statistical error. 



in a better estimate of the average distance modulus with 
redshift. This correction is based on simulations of the 
SDSS-II SN Survey, and thus is only valid for these data 
samples and selection criteria. 

We note that the differences between the plot of A/i 
in the bottom right panel of Figure [17] and that in Fig- 
ure [TO] are the reflection of more than simply an increase 
in statistics. In the public SNANA simulations of KlOa 
there existed an offset of 0.27 mag between the brightness 
of SNe la created in the simulations and those observed 
in iSullivan et al.l (|2011| ) , resulting in fainter SNe la in 
the public simulations than are used in Figure 1191 Thus 
the public simulations, while useful for our contamina- 
tion analysis, should not be used for detailed cosmologi- 
cal calculations, although we have confirmed (with newer 
simulations) that Figures [4] and [5] remain unaffected by 
this offset. 

The Malmquist bias we find is larger than has pre- 
viously reported by other surveys. This is primarily be- 
cause we are pushing the SDSS-II SN survey to its limit of 
low S/N observability; by z = 0.5 only the very brightest 
SNe are observed. We note that the ESSENCE survey 
also found a significant Malmquist bias via their simula- 
tions, which they corrected for by adjusting the prior on 
the host galaxy extinctio n (A v ) as a function of redshift 
(|Wood-Vasev et al.ll2007ft . This is an alternative method 
to adding an average redshift-dependent correction to 
each SN la distance modulus, presented here. 

4.1.1. Testing the Malmquist bias 

The above Malmquist bias correction is computed 
using a particular set of cosmological parameters. To 
determine whether this assumption biases our results, 
we have re-run our simulations with a set of widely 
varying input cosmologies, including (fl m ,fl\,w) = 
(1,0,-1), (0.5, 0.5, -1.5), and (0.2,0.8,-0.8). We find 
that the parameters describing the Malmquist bias cor- 
rection in Eqn [4] do not change beyond their quoted sta- 
tistical uncertainties. 

The value for a we have used in the simulations for 
deriving the Malmquist bias correction (a=0.11) is much 
lower than what we recover from the data (a=0.22; Sec- 
tion [63]). We examine the effect this has on our results 
by determining the Malmquist bias correction from a new 
set of simulations that uses a=0.22 in the input model. 
We find the resulting Malmquist correction is consistent 
with what we previously derived; the b parameter of 
Equation 0] is only changed by 0.14 (less than half the 
error), and the a and c parameters are unchanged. 

To demonstrate the expected effect of the Malmquist 
correction on our photometrically-classified SN la sam- 
ple, we draw ten subsamples of the same size and 
redshift distribution as our real data at random from 
photometrically-classified SNe la in our simulations. We 
derive the best-fit cosmology for these samples both with 
and without the Malmquist bias correction, again us- 
ing CosmoMC. For the uncorrected case we again find a 
best-fit value of w = — 0.87 ± 0.03, which is the weighted 
mean and uncertainty from the 10 randomly-drawn sam- 
ples. This result is significantly different from the input 
cosmology of w = — 1, demonstrating the importance of 
the Malmquist correction, the application of which pro- 
duces the expected best-fit value of w — —1.00 ± 0.03. 



17 



4.2. Host Galaxy Followup Bias 

An additional source of potential bias is due to the 
spectroscopic follow-up program of SN host galaxies. 
Our ancillary BOSS spectroscopy of likely SN la host 
galaxies (Section I2.2.1[) had an apparent (fiber) mag- 
nitude limit of rfibcr = 21.2. This should not cause 
an additional Malmquist bias, as the target selection 
is biased against fainter galaxies, not fainter SNe. It 
is well known, though, that faint galaxies preferentially 
host only luminous (high Ai) SNe la (jHamuv et al.l 
119961: IGallagher et al.ll2005D . and thus the target selec- 
tion would appear to favor detection of fainter SNe la. 
However, it has recently been shown that SNe in massive 
galaxies tend to be over-luminous for their light-curve 
shape and color (IGallagher et al.ll200l iKellv et al.l l201d: 
iSullivan et all 120101 : iLampeitl et al.l I2010H ). These ef- 
fects combine with our host-galaxy magnitude limit in a 
complicated manner that is not captured in our simula- 
tions, but which merits further study in the future. 

It is important to remember that despite these po- 
tential biases, photometric classification yields a less bi- 
ased host-galaxy sample than our spectroscopic sample. 
In Figure [30] we show the color-magnitude diagram for 
the host galaxies in our BOSS sample, with the sub- 
sample of spectroscopically-classified SNe la shown in 
blue. For each galaxy, we plot the g — r color and 
the absolute r-band model magnitude, both of which 
have been k-corrected using the standard SDSS soft- 
ware (|Blanton fc Roweisll2007D . As default, we quote all 
absolute magnitudes and the g — r color at z = 0.1. 
The host galaxies of the spectroscopic subset are, on 
average, fainter than the whole population of BOSS 
host galaxies. A Kolmogorov-Smirnov (KS) test of the 
k-corrected model r— band absolute magnitude distri- 
butions of the two galaxy samples (spectroscopically 
confirmed SN hosts and photometrically confirmed SN 
hosts) confirms the two distributions are not the same 
at a statistical significance of 99.9%. The photomet- 
ric sample includes intrinsically brighter host galaxies, 
which may be due to the increased volume sampled by 
the photometric sample as such luminous, massive galax- 
ies are rare. It could also be a product of the SN spectro- 
scopic follow-up avoiding the brightest hosts, as in these 
cases it is more difficult to separate the SN from the host 
galaxy light. However, there does not appear to be a bias 
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model colors, which is reassuring. 



5. HUBBLE DIAGRAM 



We show in Figure [21] the Hubble diagram for our 
final sample of 752 photometrically-classified SNe la, 
derived from the SDSS-II SN Survey photometry and 
SDSS-III BOSS host-galaxy spectroscopy (Section ET2|). 
In contrast to Figure I16| we have now applied our 
Malmquist bias correction to the Hubble diagram as de- 
rived in Section I4TT1 We have not corrected our sample 
for host-galaxy mass correlations, as it is beyond the 
scope of this paper (see Section l7TT|) . 

For comparison, we highlight in Figure [21] the sub- 
sample of 208 SNe la in our photometric sample that have 
spectroscopic confirmation from the SDSS-II SN Survey 
(shown in blue), and label this subsample "spec la". 
Therefore, 544 of our photometrically-classified SNe la 
are have no spectroscopic information at all, comprising 
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Fig. 20. — The host-galaxy g — r color as a function of ab- 
solute r— band model magnitude (both quantities k-corrected to 
2 = 0.1) for the host galaxies of the photometrically-classified 
SNe la (black dots) and the subsample of SNe la that have been 
spcctroscopically-confirmcd (blue cross symbols). 



72.2% of the sample. We note that only 115 of these 544 
SNe la have been previously photometrically-classified, 
using host-galaxy spect ra from the SDSS-I/II surveys 
(Sll: lHlozek et alJl2011l ). The data for all SNe la in our 
final sample can be found in Appendix [E] 

5.1. Increased scatter 

The bottom panel of Figure [21] appears to show 
an increase in the scatter of the Hubble residuals 
for the photometrically-classified SNe la (black points) 
compared to the spectroscopically-confirmed subsam- 
ple (blue points). In Figure I22| we study this appar- 
ent increased scatter by comparing the distribution of 
Hubble residuals (A/x = /i Q b s — /Avmap) in the "spec 
la" subsample to our full photometric sample, assum- 
ing the latest WMAP+BA O+g best-fit cosmological 
model (jJarosik et al.l l20TTh ■ We show these residuals in 
three redshift bins of width Az = 0.1 over the redshift 
range 0.1 < z < 0.4, which corresponds to the range 
of redshifts where these two sets of SNe la significantly 
overlap. The blue histograms show the "spec la" sub- 
sample (208 ), while the purple histograms includes the 
full 752 photometrically-classified SNe la, i.e., blue plus 
black points from Figure |2"T1 

To quantify the trends we sec in Figurc l2T1 (increased 
dispersion of the photometric SNe la with redshift), we 
fit the HR distributions with Gaussians, and report their 
full width at half maximum (FWHM) and centroids in 
Table [3] We fit Gaussians to avoid our analysis being 
adversely affected by the small but noticeable tails in 
these distributions, which are likely non-la SN contam- 
inants and are clearly offset from the Hubble diagram 
(Figure l2"Tj) at z~ 0.15. 

In Table [3] we see that the centroids of the best- 
fit Gaussians to both samples of SNe la are consistent 
with zero, showing no bias in their HR distributions 
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Fig. 21. — The H ubbl e diagram of the photometrically-classified SDSS-II SN la sample. We have corrected for the Malmquist bias as 
discussed in Section |4. II We use the best— fit values of a and /3 (see Section [6} and assumed the same M as in the simulations (A/=29.8) 
when creating this Hubble diagram. The SN intrinsic dispersion has been included in the error bars shown. Blue points show the subsamplc 
of SNe la that have been spectroscopically-confirmed as part of the SDSS-II SN Survey, while the black points only possess a photometric 
classification. The bottom panel shows the Hubble residuals of the data from the best-fit cosmology model (Section [Sj . 



as a function of redshift. However, the FWHM of the 
best-fit Gaussian does increase with redshift for the full 
photometric sample, and is additionally larger than the 
FWHM of the "spec la" subsample. The quoted errors 
in Table [3] are given by GAUSSFIT in IDL, and have been 
confirmed through bootstrap resampling of the distribu- 
tions (with replacement). We do not report the FWHM 
for the high redshift bin of the "spec la" subsample as it 
is unreliable due to small number statistics. 

We investigate why this trend appears in the full 
photometric sample, but not in the "spec la" subsam- 
ple. We plot in Figure l23l the maximum r-band S/N (at 
any epoch) for each SN la in our observed (left panel) 



and simulated (right panel) photometrically-classified 
samples. It is clear from the left panel of Figure [23] 
that SNe in the "spec la" (blue) subsample, at a given 
redshift, possess systematically higher S/N light-curves 
than photometrically-classified SNe la that weren't spec- 
troscopically observed (black). The average S/N for the 
"spec la" subsample is 27.4, whereas the SNe la with only 
photometric classification have an average S/N of 9.6. 
This is of course expected, as the SDSS-II SN spectro- 
scopic follow-up observations preferentially selected SN 
candidates that were easier to observe, naturally lead- 
ing to a bias in S/N for the spectroscopic sample. Thus 
the "spec la" subsamplc has a smaller scatter because 
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Redshift 
Bin 


Sample 
Type 


Gauss 
FWHM 


ian Fit 

Ccntroid 


Number 


0.1<z<0.2 


spec la 
photo la 


0.377±0.006 
0.413±0.008 


-0.008±0.008 
-0.012±0.013 


91 
124 


0.2<z<0.3 


spec la 
photo la 


0.366±0.010 
0.524±0.010 


0.010 ±0.012 
0.000±0.012 


80 
249 


0.3<z<0.4 


photo la 


0.610±0.016 


-0.005±0.016 


251 



TABLE 3 

Parameters of the best-fit Gaussian distributions to the data shown in Figure \TI\ 



Photometric la (753) 
Spec la (209) 
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Fig. 22. — The distribution of Hubble residuals as a function of 
redshift for the "spec la" subsample (blue) and full photometric 
sample of 752 SNe la (purple). 

it contains the brightest SNe la from the whole popula- 
tion, which are then easier to fit and thus produce tighter 
distance modulus constraints. At z > 0.3, we see in Fig- 
ure [23] the emergence of an apparent detection limit at 
S/N ~ 4—5. We have determined that this limit is 
due to the AVcolor cut (Section I3.1.5J1 . as the SALT2 
parameters are not well-determined for such noisy light 
curves, frequently returning unphysical derived parame- 
ters which we then exclude. 

6. COSMOLOGY ANALYSIS 

6.1. Fitting issues 

To allow consistent comparisons with SNLS, we use 
the same two met hods to perform cosm ological fits to 
our da ta as used in IGuv et all (|2010l ) and iSullivan et al.l 
(pcTTTh . First, we use a grid-based search technique 
when fitting simple cosmological models to just our 
photometrically-classified SN la data. Specifi cally, we 
use the simplejcos fitter (|Conlev et al.l 1201 1[ ) software 
package, which computes the \ 2 & t every point in a reg- 
ular grid (101 by 101) and converts those measurements 
to a probability via P « exp(— x 2 /2), where the propor- 
tionality is set by normalizin g over the g rid. As in the 
SNLS analysis of IGuv et al.l (|2010j ) and ISullivan et all 
we then marginalize over the SALT2 SN param- 
eters (a, {3, M) to generate confidence contours for the 
cosmological parameters of interest. 
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Fig. 23. — Maximum observed r-band S/N (at any epoch) from 
both the observed (left-hand panel) and simulated (right-hand 
panel) light-curves, as a function of redshift. The data sample is 
divided into SNe la with (blue) and without (black) spectroscopic 
confirmation. 



We also use the CosmoMC (|Lewis fc Bridle) I2002T) 
software package when fitting more complex cosmological 
models to our data (and simulations; see Section [3.1. II) in 
combination with other cosmological information. This 
package uses the Markov Chain Monte-Carlo (MCMC) 
technique to efficiently probe multi-dimensional param- 
eter space, allowing one to quickly investigate a large 
number of different regions in the parameter space. We 
have slightly modified CosmoMC to allow for the simul- 
taneous fitting of both the cosmological parameters and 
the SALT2 SN parameters a and (3, which define the 
standardization of SNe la. We include in the distance 
modulus calculation the redshift-dependent Mahnquist 
bias correction (Eqn [4| and the full SALT2 light-curve 
parameter covariance matrix. Finally, at each point in 
the MCMC chain, we determine M (absolute magnitude 
at peak for SNe la) for the value of H at that step 
in the chain. This approach is the sam e as analyti- 
cally marginalising over M, as outlined in iBridle et al.l 
( 20011) . and is a sim ilar methodology as used by SNLS 
([Sullivan et al.l 120 111 ) . This method is used in all our 
cosmological fits. 

We execute our modified CosmoMC code using six 
chains in parallel to facilitate quick coverage of the 
large, multi-dimensional parameter space. Each chain is 
started at a random location within our defined paramc- 
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ter space and typically converges after 50,000 to 100,000 
steps. We assume that the MCMC has converged when 
R— 1 < 0.1, where R is the Gelman and Rubin statistic, 
i.e., when the variance w ithin the chains is equal to the 
variance between chains (jBrooks fc Gelmanlll998| ). 

We provide in Table 4 a summary of the different 
combinations of priors we assume during our different 
cosmological fits discussed below. Throughout all our 
analyses we assume flat distributions for the priors on our 
SN parameters, i.e., a = (0.01,0.5) and j3 = (1.0,5.0). 
We also assume flat distributions for the priors on fid m 
(density of dark matter) and fib (density of baryonic mat- 
ter) when these cosmological parameters are allowed to 
vary. 

Finally, we set the value of the intrinsic dispersion of 
SNe la (cri n t) to C7i n t = 0.12 mag, which is then added in 
quadrature to all our SN errors. Although this results in 
a reduced \ 2 close to unity for all our cosmological fits, 
the best-fit value (i.e., delivers the reduced \ 2 closest to 
one) for all our data is 0.16 mag. However, our simu- 
lations show that this larger value of <Ji n t is consistent 
with our small level of non-la SN contamination^], i.e., 
we input <7i n t = 0.12 mag in our simulation, but measure 
Pint = 0-16 mag for a photometrically-classified sample 
like ours. As ai n t is used to explain unknown residual 
scatter in the SNe la population, we feel it is appropriate 
to remove any extra scatter caused by non-la SN contam- 
ination from our 'measured' value of a- m t. We stress that 
this statement is not in contradiction with Section 15.11 
where we state that the larger scatter we observe in our 
sample could be caused by lower S/N light curves, as such 
lower-quality data would increase both the scatter in the 
population and the observed errors bars on our SN dis- 
tance moduli. However, contamination would increase 
the overall scatter in the SN Hubble residuals without 
necessarily increasing the distance uncertainties. It is 
worth noting that we also ran our cosmology fits with 
the (Tint =0.16 mag and found no significant difference 
to our cosmological results. 

We note that our assumed value of (Tint = 0.12 
mag is still higher than that fo und in ILampeitl et al.1 
(|2010al cri„t =0.088 mag) and in iKessler et al.1 (|2009aL 
cr int =0.08 mag, SALT2 and SDSS-II SNe). These val- 
ues were computed for the first-year SDSS SN Sur- 
vey spectroscopic sample; limiting ourselves to the 
spectroscopically-confirmed subsample of our full pho- 
tometric sample (Section 5), we still fin d a larger in- 
trinsic scatter, c^nt =0.11 mag. However. iConlev et al.l 
(12011ft found o-int = 0.10 mag for the SDSS spectro- 
scopically confirmed data, in agreement with our value. 
The apparent discrepa ncy may be caused by K09 and 
ILampeitl et al.l (l2010a[) using an old er version of SALT2 
(|Guv et al.1 12007ft . as IConlev et all (120111) also us es the 
same newer version of SALT2 (iGuv et al.l 12010ft used 
here. As our simulations have shown that our selection 
criteria do not result in a bias between the input and 
best-fit value of a- ln t for a pure (i.e., spectroscopic) sam- 
ple of SNe la, we have confidence in and continue to use 

6 We find that the intrinsic dispersion of our full sample drops 
to (Tint = 0.1, and our cosmological results remain consistent, if we 
simply remove the 25 SNe la that are located > 3<r away from the 
best-fit cosmological model. We do not recommend such a "sigma- 
clipping" technique when using SNe la to test cosmological models, 
but this test does illustrate the sensitivity of ai n t to such outliers. 
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TABLE 4 

Priors imposed on the fitted cosmological parameters in 
four different combinations (sets). the dash symbol in the 
table represents where we do not need to set priors for 
parameters, as they are constrained by a combination of 

OTHER PRIORS (E.G., f2 m IS RESTRICTED BY THE PRIORS ON Q dm 
AND f! b ). 

throughout our value of cri nt = 0.12 mag. We stress that 
the assumed value of <Ti nt , within the range discussed 
above, has little effect on the cosmological fits presented 
in this paper as the fits also have the freedom to adjust 
the values of other SN parameters, such as a and /3. 

6.2. Our constraints on ACDM 

We begin by studying the cosmological constraints 
obtained using only our SN la sample, before combining 
with other data. For completeness, we provide a sum- 
mary of all our cosmological fits, including the different 
combinations of data-sets and priors, in Table [5] 

We begin by fitting the ACDM cosmological model 
(w = —1) to our photometrically-classified SNe la us- 
ing simple -cos fitter and prior Set I in Table [4] Under 
the assumption of flatness we obtain a best-fit value of 
fi m = 0.24^q'q5 (statistical errors only). When we relax 
the prior on flatness (prior Set II), we obtain the (grey) 
confidence contours for fi m and fl\ in Figure [24] For 
comparison, we show similar constraints on these cos- 
mological par ameters us i ng th e three-year SNLS data 
(SNLS3) from lGuv et "all (f2010h . which only includes 242 
spectroscopically-classificd SNe la from SNLS. 

Figure [M] also demonstrates that our 
photometrically-classified SN la sample alone is able 
to detect an accelerating universe (i.e., fiA > fi,„/2). 
Integrating over the whole parameter space in Figure [21] 
(see table [4] for parameter ranges), we compute the 
probability of an accelerating universe given our data to 
be 99.96% (statistical uncertainties only). 

6.3. Our constraints on wCDM 

We next fit for a flat, wCDM cosmological model 
using CosmoMC and the prior Set III in Table [4] Q. 
We fit this model to our sample of photometrically- 
classified SNe la and the Ho measurement of 73.8 ± 
2.4 km s _1 Mpc~ x from the recent "Supern ovae and Hp 
for the Equ ation of State" sample (SHOES; iRiess et al.l 
120091 12011ft . This Gaussian prior on Hq does not impact 
the cosmological constraints due to our marginalization 
over the absolute magnitude of SNe la, but rather en- 
sures that cosmoMC, which assumes a minimum age to 
the Universe, performs in a well-behaved manner. We 

7 All other cosmological parameters are left at their default val- 
ues (i.e., re-ionization optical depth, the primordial super-horizon 
power in the curvature perturbation on 0.05 Mpc scales, and the 
scalar spectral index). 
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Fig. 24. — The 68% and 95% confidence contours of f2 m versus Q\ for a ACDM model using only our photometrically-classified SNc, 
with prior Set II (Table 5), and allowing curv ature to vary. On ly statistical errors in the contours are shown. The blue dashed contours 
show the comparable SNLS3 constraints taken Guy ct al. (2010). 

note that within cosmoMC, f2 m is a parameter comprised 
of two components: fib and Odm. Since the ratio of 
these two components is not constrained by SNe la alone, 
we fix n h to the WMAP7 value of 0.0458 (jJarosik et al.1 
120 111 ), reducing the number of parameters to be con- 
strained by the SNe data to two: fidm and w. 

We find the best-fit values for this model to be 
w = -0M±°H and fi m = 0.2lt° \l (statistical errors 
only), and show in Figurel25lthe joint 68% and 95% confi- 
dence intervals for this data. We also show the equivalent 
contours for the three-year sample of spectroscopically- 
confirmed SNe la from the SDSS-II SN Survey. For this 
comparison, we include all 306 confirmed SNe la (regard- 
less of whether they are part of our photometric sample 



or not) at a redshift of z < 0.3, below which any selection 
bias in the spectroscopic sample should be minimized. 
This comparison demonstrates that using photometric 
instead of spectroscopic classification, which increases 
the size of the sample by a factor of 2.5 and extends 
the redshift range, results in a reduction of the area of 
the confidence contours by a factor of 1.6. We stress that 
this simplistic comparison does not constitute a detailed 
analysis of the full 3-year spectroscopic sample of the 
SDSS-II SN Survey, which will be presented elsewhere. 
We also find that limiting our photometric SNe la sam- 
ple to z < 0.3 gives consistent cosmological constraints 
with the samples plotted in Figure 1251 but with slightly 
larger uncertainties. 
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TABLE 5 

Summary of the cosmological fits presented in Section [6] 



During this analysis, we simultaneously solve for the 
best-fit values of the SALT2 SN parameters. In our 
fit of the wCDM model, we find a 



0.22™ and 



p = 3.12_ 012 . Our fitted value for (3 is in agreement 
with previous analyses of the SPSS data (iLampeitl et all 
l27)10al IMarriner et all I27LTT1 : [Conlev et al.l 120 lift . How- 
ever our value of a appears higher than previous anal- 
yses; [Lampci tl et al.l (|2010aD found a=0.16±0.03, and 
IMarriner etldl (|201lD found a=0.13lt^t One poten- 
tial explanation of this difference could be the non-la 
SNe contamination in our photometric sample, but this 
has been tested in our simulations: wc recover the in- 
put a and (3 both with and without the expected level 
of contamination detailed in Section [3J A more plausible 
explanation appears to be the higher average S/N of the 
spectroscopically-confirmcd SNe la; using only the sub- 
set of spectroscopically-confirmcd SNe la in our photo- 
metric sample, we find a=0.16±0.02, in agreement with 
previous results. 

We explored whether HR outliers (which we do not 
automatically clip) might be affecting the derived value 
of a. We subtracted a * Xi from the HR of each SN 
in our sample and plotted these as a function of X\, 
finding four clear outliers. Removing these objects from 
our sample and refitting the cosmology, wc found the 
resulting cosmological contours to be unchanged but the 
derived SALT2 parameters to both be lower (a=0.18, 
/3=2.79). As this a is still larger than what is found in 
previous studies, the outliers cannot be solely responsible 
for this discrepancy, and regardless does not result in 
biased cosmological constraints. 

6.4. Constraints from combining data-sets and 
comparison with SNLS 

Finally, we determine cosmology constraints with 
our photometrically-classified sample of 752 SNe la com- 
bined with cosmological information from the power 
spectr um of Luminous R ed Galaxies (LRGs) in the SDSS 
DR7 (IReid et al.l [20loh the full WMAP7 CMB power 
spectrum (|Larson et al.ll2011f) . and the SHOES H mea- 
surement. The SHOES Hq measurement is partially de- 
termined using nearby SNe la measurements, and thus to 
be fully consistent we would have to consider the covari- 
ance between this value of Hq and our SNe la measure- 
ments. However, as we are assuming no prior information 
on M in our treatment of intrinsic SN parameters, these 
measurements can be considered independent. Further- 
more, the uncertainty in M is a subdo minant systemati c 
uncertainty to the derived value of Hq (|Riess et al .1120 111 ) . 

We fit this combination of data using CosmoMC 
for a non-fiat wCDM cosmology, using the priors listed 
as Set IV in Table SJ With the addition of these ex- 
ternal data-sets, we can now relax our priors on the re- 
ionization optical depth (t=[0.00, 0.50]), the primordial 
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Fig. 25. — The 68% and 95% confidence contours of w and f2 m in 
a flat uiCDM model (assuming prior Set III in Table 5) for a combi- 
nation of our photometrically-classified sample of 752 SNe la and 
the SHOES measurement of Ho (only the statistical errors shown, 
grey ) . The blue contours show the equivalent constraints but using 
the 3-year sample of spectroscopically-confirmcd SNe la from the 
SDSS-II SN Survey (confined to z<0.3). 



super-horizon power in the curvature perturbation on 
0.05 Mpc -1 scales (log A= [0,30]), and the scalar spectral 
index (ro s =[0,1.5]). 

We find the best-fit value for the equation-of-state 

+0.10 
-0.10' 



of dark energy using these data is 



-0.96 



with n m = 0.291^2 and fi k =0.00l£^ (statistical errors 
only). Wc also find a best-fit value of H =67 '.97±%™ 
(stat) km s _1 Mpc -1 . These cosmological constraints 
are summarized in Table [SJ 

Figure [21)1 shows our joint confidence contours for w 
and f2 m in compari s on wi th similar SNLS3 constraints 
from iSullivan et all pOlTh using the same combination 
of external data-sets (CMB, LRGs, and the H SHOES). 
SNLS3 only uses spectroscopically-classified SNe la, col- 
lected from the following SN data-sets: 242 SNe la from 
SNLS; 123 low-redshift SN e la from the literature (pri- 
marily iHamuv et al.l 119961 iRiess et al.l Il999t [Jha et al.l 
120061: iHicken et al.ll2009t iContreras et alj|200£ ): 14 high 
redshift SNe la from HST (|Riess et~alTl2007T) : and 93 
SNe la over the redshift ran ge 0.06 < z ^ 0-4 fro m the 
first-year SDSS SN Survey (jKessler et aLll2009al) . The 
inclusion of the SDSS SN data means there are some 
SNe la in common between these analyses. There ap- 
pears to be good agreement in these two sets of con- 
straints, although there is a small offset in the best fit 
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Fig. 26. — The 68% and 95% confidence contours of w and 
Q m in a uiCDM mod el from the spectroscopic SNLS3 sam- 
ple HSullivan et aI.ll20Tll . blue) and our photometrically-classified 
SDSS-II SN la sample (grey) using prior set IV in Table 5. Both 
sets of contours include external data from the CMB, LRGs, and 
the SHOES Ho measurement. Contours represent statistical uncer- 
tainties only. 

values for f2 m between the two analyses: Af2 m =0.018, 
significant at less than la. The best-fit values of w differ 
by Aui=0.080 between the two analyses, which is again 
consistent within the quoted uncertainties. One possible 
explanation for these small differences is the fact that we 
have not corrected for the correlation with host-galaxy 
stellar mass, as discussed in Section 7. Overall this com- 
parison is reassuring, considering the lower redshift lever- 
age of the SDSS-II SN sample (z < 0.55) and the lack 
of spectroscopic confirmation used herein. These results 
demonstrate the potential for photometrically-classified 
SN la samples to be used to improve cosmological con- 
straints in the future. 

7. DISCUSSION 

7.1. Systematic Uncertainties in Photometric 
Classification 

The goal of this paper has been to illustrate the 
power of photometrically-classified SNe la in deliver- 
ing competitive cosmological results when compared to 
spectroscopically-confirmcd samples. As such, we have 
paid particular attention to the sources of uncertainty 
unique to our methodology of SN classification, and to 
the SDSS-II SN Survey in general, i.e., contamination 
from non-la SNe and Malmquist bias. We have not un- 
dertaken an exha ustive study of system atic uncertainties 
as exe mplified bv lKessler et al.l (|2009aD and lConlev et al.l 
(|20ll . 

In Section 21 we studied the effect of non-la SN con- 
tamination of photometric samples on cosmological con- 
straints. The low predicted contamination rate (3.9%) 
in our sample has an insignificant effect on the best-fit 
cosmological parameters in simulations; compared to a 
completely pure sample of SNe la we find a bias in the 
cquation-of-state of only w ~ 0.007 (using prior Set III 



in Table 0}. 

We have also studied the systematic offset associ- 
ated with the Malmquist bias effect (Section 14. ip . This 
correction, as a function of redshift, is included in all 
our cosmological fits. We also investigated allowing the 
three parameters in the Malmquist bias exponential pa- 
rameterization (Eqn 4) to vary in the cosmological fit 
within the error bars; the cosmological constraints re- 
mained identical. Correlations between the three param- 
eters in the Malmquist bias fit have not been investigated 
in this paper. 

We do not attempt to correct for the known cor- 
relation between the SN la Hubble residuals and the 
properties of their host galaxies; SNe la in massive 
galaxies are over-luminous even after light-curve cor- 
rections. Our sample of host galaxies span a range in 
absolute magnitude of approximately —24 < M®' 1 < 
— 18 (Figure 20), which corresponds to a stellar mass 
range of ~ I0 9 -I0 11 M^, or a p redicted difference of 
Afj. ~ 0.1 (jLampeitl et al.ll2010bf ). We additionally ex- 
pect that the magnitude limit imposed on our host- 
galaxy selection could cause a bias by preferentially 
selecting more massive (more luminous) host galaxies 
at higher redshift, thus preferentially selecting over- 
luminous SNe la. Unfortunately, the underlying phys- 
ical mecha nism that drives this correlation remain s un- 
clear (e.g. JD'Andrea et al.llMfllGupta et al.ll20Ifl ). and 
as such the mass-HR relationship could also be subject 
to a redshift dependence. Therefore, we do not correct 
for this effect here, but note that it will be essential for 
future studies. We note that introducing a correction 
term for host-galaxy properties is no more difficult in 
photometrically-classified SNe la samples than in those 
obtained from spectroscopic follow up, so it is a common 
problem for all future SN surveys. 

A systematic uncertainty that is unique to 
photometrically-classified SN la samples like ours, based 
on host galaxy spectroscopy, are errors in associat- 
ing the SNe with the correct host galaxy. As de- 
scribed in Appendix [C] we found only one mismatch 
between the SDSS-II SN spectroscopic subsample and 
our photometrically-classified SNe la, one mismatch be- 
tween the SDSS-II host galaxy spectroscopy and BOSS 
host galaxy spectroscopy, and giving an error rate of 
only 0.6%. We have removed this object, as well as an- 
other likely matching error discussed in Section[X2J from 
our photometric sample. We cannot rule out the possi- 
bility that other SN-galaxy pairs have been incorrectly 
matched, although based on the rate found in the spec- 
troscopic subsample we expect the number to be low and 
thus have a negligible effect on our cosmological fits. 

Finally, our estimation of systematic effects due to 
Malmquist bias and contamination relies on the assump- 
tion that our simulations accurately represents the final 
sample after selection cuts. However, some implicit as- 
sumptions in the simulation may not be valid. For ex- 
ample, we assume that all candidates are SN la or non- 
la SNe, even though some photometrically-classified can- 
didates may in fact be another type of transient (e.g., 
AGN). We also assume that the 41 non-la templates 
reflect a complete sample of non-la SNe, and that the 
non-la properties are redshift independent. These as- 
sumptions may be inadequate, as discussed further in 
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Section 17^1 and iKessler et al.l (|2010aD . 

The magnitude of these systematic biases on our 
results are either small compared to our statistical un- 
certainties or unknown, and as such wc do not include 
them in the error budget of our derived cosmological pa- 
rameters in Table [5j 

7.2. Future Improvements to Photometric Classification 

We have provided in this paper a procedure for 
photometrically-classifying SNc. Of course, other meth- 
ods also exist, with their own relative advantages and 
disadvantages. For example, rather than applying limits 
on light-curve properties in 2-D parameter space (e.g., 
Figure [6]), one could apply a nearest- neighbor algorithm 
to look for clustering in higher dimensional parameter 
space. Alternatively, we note that the analysis presented 
here excludes from our sample all transients that have 
Pi a computed to be below a hard threshold (see Sec- 
tion [3113). One could choose instead to retain all SN 
candidates in their cosmology analysis, weighting each 
candidate by its Bayesian probability (Pi a ) of being a 
SN la. This approach avoids the uncertainty of choos- 
ing the optimal P\ a threshold to obtain a "clean" sample 
of SNe la, and prevents the removal of actual SNe la 
and the information that they provide. These methods 

- the nearest-neighbor algorithm an d the fully Bayesian 
method using the BEAMS algorithm (jHlozek et al.ll2011ft 

- are currently being investigated using the SDSS-II SN 
candidates (and their BOSS host-galaxy redshifts). 

Although we have optimized the selection cuts im- 
plemented in this paper, utilizing higher-order criteria 
could further improve the accuracy and efficiency of pho- 
tometric classification. For example, the ellipsoidal X\- 
color cut could be derived with an additional rotation 
parameter that accounts for the correlations between the 
two SALT2 light-curve parameters. 

There is also room for improvement in the simu- 
lations we use to describe our observations and measure 
the efficacy of our selection criteria. The SN la templates 
included in the simulations are missing a number of sub- 
classes, such as 02cx-like SNe, which have sim ilar light 
curves to but are fainter than normal SNe la (|Li et al.1 
12001 IPhillips et afll2007l ; IMcClelland e t al.l l2010fl: 06bt- 
like SNe, which are particular ly problematic ( Foley et all 
l2010al iStritzinger et al.ll2011f); and super-Chandras ekhar 
mass SNe la ( Howell et al 

simulations also include only a limited number of non- 
la SN templates, a deficiency which is sure to be im- 
proved in the future as the quantity of observations of 
these objects begins to reflect their diversity. 

At present, the only transient objects included in 
these simulations are SNe (la, II, and Ibc) . In the future 
we would like to include other transient objects that are 
known contaminants in SN surveys, such as AGN. There- 
fore, when we quote an estimated contamination of 3.9%, 
this could be underestimated as we have only included 
contamination from non-la SNe. Furthermore, errors on 
assigning the correct SN host redshift were not modeled 
in our simulations. This is a rather complicated effect, as 
it depends on both the distribution of SNe within galax- 
ies and the luminosity function of galaxies (those with 
and without SNe) as a function of redshift. While we be- 
lieve we have removed possible mis-identified host galax- 
ies from our data (Appendix [C]) , it would be interesting 



to model whether the expected number of such misiden- 
tifications is consistent with our findings. 

In this work we have assumed a constant value 
for the intrinsic di spersion of our sample. However, 
IKessler et al.1 ()2012T ) have recently shown that the intrin- 
sic dispersion of SDSS and SNLS SN data may be better 
described wavelength-dependent function. There- 
fore, future analyses will want to include a more com- 
plex model than is currently standard for the intrinsic 
dispersion to improve their cosmological constraints. 

We note that the rela tive rates as a funct ion of red- 
shift of SNe la (base d on iDildav et~al~l l2008| ) and non- 
la SNe ([B azin et all l2009t) assumed in the simulations 
still have associated uncertainties. It is certain that these 
constraints will be improved in the near future by the 
next generation of large, deep SN surveys. 

Finally, in this paper we have used a simple redshift - 
dependent correction for the Malmquist bias. In the 
future, higher-order Malmquist bias corrections could 
be investigated, such as stretch-dependent or color- 
dependent corrections. The latter of these may be im- 
portant, as Figure fT5l shows there is a clear bias in the re- 
covered color distribution with redshift. The ESSENCE 
survey used a color-dependent Malmquist bias correc- 
tion, adjusting the prior on t he host galaxy extinctio n 
(A v ) as a function of redshift dWood-Vasev et aLll2007l) . 
This method would be interesting to explore with our 
photometrically-classified SNe la. 

7.3. Prospects for Future SN Surveys 

The DES Supernova Survey (|Bernstein et al.ll201lD 
should start by the end of 2012 and run for at least 
five seasons. It is expected to measure high-quality light 
curves for ~ 4000 SNe la out to a redshift of z w 1.2, 
for which real-time spectroscopy of every SN la can- 
didate, as done previously, will be impractical. The 
SNLS, SDSS and ESSENCE surveys combined used over 
a year of telescope time (4 and 8-m clas s) to spectroscop- 
ically confirm fewer than 1000 SNe la (|Folev et al.ll2009t 
IHowell fc Legacy Survev l I2009I ). 

Therefore, DES will need to use photometric classi- 
fications to reduce the burden of real-time spectroscopic 
confirmation and use allocated spectroscopic resources 
wisely, e.g., targeting hostless SNe for spectroscopic con- 
firmation or building up training samples of SNe for pho- 
tometric classifiers. Our work suggests the need for ob- 
taining spectroscopy of SN host galaxies which, as was 
the case of SDSS-II, can be done over a longer period 
of time and can be c oordinated with other science goals 
(|Lidman et al.ll2012ft . 

Photometric classification is well suited for obtain- 
ing large, uniformly-selected samples of SNe la. In 
light of the systematics discussed in Section 17.11 and Ap- 
pendixini there is a great scientific benefit in having such 
samples that can be subdivided and analysed a number 
of different ways to determine the magnitude of such sys- 
tematic effects. We describe below two examples of how 
these samples will be useful. 

SN lensing is the increase in observed flux from a SN 
due to lensing by the structure the light passes through 
on its journey through the universe. iClarkson et al.l 
(|2011f ) has discussed how the number of galaxies we ob- 
serve along the line of sight to a SN should be correlated 
with the over-luminosity of the SN (voids play a role in 



25 



this determination as well). Thus grouping large num- 
bers of SNe la by the amount of foreground structure 
is vital for measuring this effect; such a program is cur- 
rently being investigated with this sample (Smith et al. 
in prep). 

The correlation between SN la host-galaxy prop- 
erties (mass, star-formation rate, metallicity) and 
SN la HRs after light-curve and color corrections 
are made is an importa n t disc o very for SN l a cos- 
mology (iGallagher et all 120081: iHicken et all 120091: 
Kellv et al.l 120101: iSullivan et all 120101: ILampeitl et al.l 
2010bt IGupta et all 120111: ID'Andrea et all 1201117 7 
Using spectr o scopica lly-confirmed SDSS-II SNe, 
ILampeitl etldl (|2010bD showed that SNe la are -0.10 
mag intrinsically over-luminous in passive hosts, and 
also found that parameterization using stellar mass gives 
an improvement of — 4er on the cosmological constraints. 
This effect, however, is most likely no t driven by the 
host m ass; using low-z SDSS-II SNe la, ID'Andrea et~aH 
(|2011| ) found intrinsically over-luminous SNe occur in 
high-mctallicity gala xies at >3g. The r ecent SNLS3 
cosmology analysis (jSullivan et al.l 120 111) is the first 
major SN study to include host-galaxy corrections in a 
full cosmology analysis. 

Our SDSS-II photometric SN la sample is much 
larger than any previously analysed sample, and the 
improved statistics and reduced bias may lead to 
an improved understanding of this effect. Cor- 
relations with photometrically-derived stellar mass 
and spcctroscopically-derived metallicitics and star- 
formation rates are typically degenerate, but a large sam- 
ple will allow one to hold these parameters constant while 
allowing only one to vary. 

8. CONCLUSIONS 

In this paper we use the full three-year photometry 
from the SDSS-II SN Survey, together with BOSS spec- 
troscopy of the host galaxies of transients, to create a 
photometrically-classified sample of SNe la to be used 
for cosmology. Our main results are: 

• We have created a homogeneous sample of 752 
photometrically-classified SNe la; the largest col- 
lection of SNe la ever selected from a single photo- 
metric survey. Our sample spans a redshift range of 
0.05 < z < 0.55 and contains 544 newly classified 
SNe la; 208 SNe la in our sample had been previ- 
ously spectroscopically confirmed and another 115 
had been photometrically-classified using SDSS-II 
host galaxy spectra. Based on SNANA simula- 
tions, we estimate that this sample is 70.8% effi- 
cient at detecting SNe la, with a contamination of 
only 3.9% from core-collapse SNe. We demonstrate 
that this level of contamination is negligible when 
estimating constraints on cosmological parameters 
using this sample of SNe la. 

• Malmquist bias and SALT2-rclated effects are the 
largest systematic selection uncertainties in our 
photometrically-classified SNe la sample. We esti- 
mate the combined size of these effects using exten- 
sive SNANA simulations, which show that they can 
be > 0.1 magnitudes at the high-redshift limit of 
our sample (z = 0.5; Figure [T9|) . We use our simu- 
lations to correct for these biases and show we can 



then recover the correct cosmological model input 
into the simulations. 

• We show that the "spec la" subsample of 208 
SNe la is potentially biased, as they possess 
higher S/N light curves than the majority of our 
photometrically-classified SNe la ( Figure 123]) . Fur- 
thermore, there is evidence in Figures [14] and [15] 
that the "spec la" subsample is biased to brighter 
(higher stretch) and slightly bluer (lower color val- 
ues) SNe la than the whole population. The 
weighted means of the SALT2 parameters for the 
"spec la" sample are X\ = 0.033 ± 0.015 and c = 
-0.021 ± 0.002, compared to X x = -0.017 ± 0.013 
and c = —0.018 ± 0.002 for the sample as a whole. 
It is additionally clear from Figure [50] that the 
host galaxies of the "spec la" subsample are biased 
towards fainter (in absolute magnitude) galaxies 
compared to the whole population of host galax- 
ies. These biases may be related, but we have not 
investigated this possibly herein. 

• We present the corrected Hubble diagram for our 
photometrically-classified SN la sample in Fig- 
ure 1191 The extra scatter seen in this diagram at 
high redshifts is likely caused by our sample in- 
cluding SNe la with lower S/N than are found in 
the spectroscopically-confirmed SNe la subsample 
(based on our simulations). We then fit this Hubble 
diagram with a straightforward ACDM (w = —1) 



cosmological model to obtain Q m = 0.24" 
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tistical errors only) for a flat Universe. If we re- 
lax the constraint on flatness we obtain the con- 
straints on f2 m and £l\ shown in Figure 1241 where 
we have detected an accelerating Universe at the 
99.96% confidence level. This figure also shows 
that our statistical constraints on these impor- 
tant cosmological parameters are comparable to 
th e recent SNLS th ree-year constraints published 
in IGuv et al.1 (l20l"ol) . 

• In Figures [25] and [26] we show our constraints on 
the equation-of-state of dark energy (w) for our 
photometrically-classified sample on its own (with 
only Ho data) and when combined with other cos- 
mological information (CMB, LRGs, Hq). We find 
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-0.951^1 and w 
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respectively 



(statistical errors only), which are consistent with 
both the SDSS-II and SNLS three-year spectro- 
scopic samples. These cosmological analyses il- 
lustrate that our photometrically-classified sample 
can deliver competitive constraints even though 
it lacks extensive SN spectroscopic follow-up and 
probes a smaller range of redshifts (compared to 
the SNLS3). 

Creating a photometric SN la sample is a fundamen- 
tally different task from creating a spectroscopic sample. 
There are various equally-valid approaches that can be 
taken into account in designing a SN la photometric- 
classification algorithm depending on the purpose to be 
fulfilled, each resulting in a sample of different size and 
composition. As our focus is on obtaining useful cos- 
mological constraints, this required prioritizing a highly 



26 



pure sample (few non-la SNc included) over a highly effi- 
cient one (few SNe la excluded). Had our intention been 
to study, for example, SN rates, then the selection pres- 
sures driving our sample construction would have been 
different. A focus on cosmology and thus purity will ne- 
cessitate the sort of strict cuts on photometric properties 
that we apply. 

Our classification is based on the technique of Sll, 
which computes the Bayesian probability of each SN sub- 
type based on the fit of the data to templates and models 
( Appendix [A"]) . However, on its own this method does not 
perform nearly as efficiently as it does when a prior on 
the redshift exists for each SN candidate (Olmstead et 
al. in prep.). We undertook an ancillary program using 
BOSS to obtain host-galaxy redshifts for a large number 
of SN candidates; this forms a crucial part of our paper. 
We have argued that obtaining the resources to spectro- 
scopically observe the majority of SNe la in future sur- 
veys is not feasible. This is because SN spectroscopy is a 
highly time sensitive (observing windows of a few weeks) 
and scattered (low density per solid angle per unit time) 
undertaking. For this reason, it is also subject to signif- 
icant selection biases. Host-galaxy spectroscopy allows 
the observer to obtain the redshift of each object much 
more efficiently: multi-object spectrographs can be used 
to sample the higher spatial density of targets at sched- 
uled dates long after the SNe have faded away. For this 
reason we believe host spectroscopy will remain a vital 
component of SN surveys in the future. 

It is not desirable for future surveys to abandon real- 
time spectroscopy completely; it will remain necessary to 
identify subtypes of SNe la, train classifiers, and study 
detailed properties. It is even possible for spectroscopic 
samples with sizes and redshifts ranges exceeding that of 
SNLS and SDSS-II to be created. However, this method 
could not possibly achieve the volume of SNe la iden- 
tification possible with photometric classification. And 
although statistical uncertainties today are quite small, 
sample size is very important for understanding system- 
atic uncertainties ( Appcndix[DJ . Future surveys will have 
to understand SN lensing, host-galaxy correlations, in- 
trinsic color, evolution, and other effects that have the 
potential to bias cosmological constraints. Large sam- 
ples that allow a complicated parameter space to be 
explored will thus be necessary. For all these reasons, 
photometric-classification as an underpinning of SN cos- 
mology is here to stay. 
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A. PHOTOMETRIC SUPERNOVA IDENTIFICATION 

In this paper we make use of the Photometric Su- 
perNova IDentification (PSNID; Sll) software to obtain 
a typing of SN candidates using only the photometry of 
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each object. PSNID has been shown to perform better 
than other photometric-classification algorithms, scoring 
the highest Figure-of-Merit in the recent classification 
challenge of KlOa. This methodology has been applied 
to the full 3-year data release of the SDSS-II SN Survey 
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and will be presented in full in Sako et al. (in prep.). 

We refer the reader to Sll for a full description of 
PSNID, but describe here some of the key features of 
the algorithm which are relevant to this paper. PSNID 
uses the observed multi-color light curve to calculate a 
reduced x 2 fit to a grid of SN la light-curve models and 
non-la SN templates, assigning probabilities for each SN 
using the Bayesian evidence criteria. PSNID returns the 
SN la Bayesian evidence (Ei a ), given by, 
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where the redshift prior P(z) for an externally con- 
strained redshift z ox t is given by, 
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and <7z is the error on the observed (external) redshift. 

The E\ a is computed by marginalizing the prod- 
uct of the likelihood function and the prior probabilities 
over the model parameter space. Five parameters are 
included in the model: redshift (z), host-galaxy extinc- 
tion (A v ), time of maximum light (T max ), the amount 
by which a SN la declined in the B-band du ring the 
first fifteen days after maximum light (Amis, IPhillipsI 
I1993D . and the distance modulus (//). Flat priors are 
assumed on Ay, T mayi , and Ami s. Galactic ex t inctio n 
is corrected for using A\r f rom ISchleeel et all (fl998l) . 
assuming the iCardelli et al.l (|1989l ) extinction law with 
Rv = A V /E(B - V) = 3.1. The color law for SN 
host gala xies, however, is assu med to be steeper, with 
Rv = 2.2 (|Kessler et al.ll2009aD . 

The non-la SNe Bayesian evidences (-Eibc,n) are 
given by, 
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This quantity is the summation over a variety of Type 
Ibc and Type II templates, as given in Table 1 of Sll. 
PSNID returns a Bayesian probability for each of the 
three SN types considered herein, given by 
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where, by definition, Pi a + Pi^ c + P\\ = 1. 

We use in this paper a version of PSNID that has 
incorporated several improvements over that presented 
in Sll. The primary difference is in the light-curve tem- 
plate uncertainties (Equations 6 & 7 from Sll); these 
have been changed to better match the observed distri- 
bution of errors in spcctroscopically-confirmcd SDSS-II 
SNe. Specifically, the magnitude errors Sm on the SNe la 
gri model light curves are now 
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0.06 + 0.04 x (|*|/20) |t|<20days, 
0.10 + 0.18 x ((|t| - 20)/60) |t| > 20 days, 

(A5) 



while the new magnitude errors for the non-la SN light- 
curve templates are 



Sm cc = 0.06 + 0.08 x (|i|/60), 



(A6) 



where t is the rest-frame epoch in days from B-band 
maximum. These changes require us to re-determine 
goodness-of-fit thresholds used in Sll on PSNID, which 
we describe in Section f3. 1.41 

In Section 12.2.11 we obtain a preliminary classifica- 
tion of our SDSS-II SN candidates by running PSNID on 
the ugriz light-curves while placing a flat prior on the 
redshift. This classification shapes our subsequent target 
list for BOSS host-galaxy spectroscopy, and is the only 
occasion in this paper where we use a flat redshift prior. 
For both our simulated SNe (Section 13. 1.4[) and our SN 
candidates that have host-galaxy spectra from BOSS we 
use the spectroscopic redshift as our prior ( Section 13. 2[) . 
This additional prior on the redshift helps compress the 
parameter space being investigated and break degenera- 
cies between non-la SNe at low-z, which could appear 
similar to higher redshift SNe la. As in Sll, we assume 
flat priors for all other PSNID model parameters. 

B. SPECTRAL ADAPTIVE LIGHTCURVE TEMPLATE 

We determine the distance moduli for our SN candi- 
dates (both data and simulated) u sing the SALT2 light- 
curve fitting software (version 2.2: IGuv et alll2010h . For 
each SN light-curve the SALT2 fitting program deter- 
mines the best-fit value of three parameters (Xq , X\ , c) , 
which describe the observed lu minosity offset, stretch , 
and color of the SN, respectively (|Guy et al.ll2007L l2010f l . 
These fitted values are then used to "standardise" the 
light curve, as the distance modulus to each SN is calcu- 
lated using, 
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where to*b is the _B-band peak apparent magnitude and 
is defined as — 2.5 log 10 (xo). Parameters a, (3 and M (ab- 
solute B-band magnitude at peak) are constants that 
can either be derived for the whole sample simultane- 
ously with the best-fit cosmology, or can be constrained 
from other data. In our cosmology analysis presented in 
Section [6] we allow a and (3 to float within priors and an- 
alytically marginalize over M (which is degenerate with 
Ho)- 

We fit (in flux) SALT2 to the SDSS SMP light-curve 
data in the griz passbands. Although the SALT2 tem- 
plate does not extend to the rest-frame z-band, we in- 
clude z-band data in our fits as it should help constrain 
the model at higher redshifts. However, the z-band has 
low throughput in the SDSS, and we obtain consistent 
light-curve fits with and without including this infor- 
mation. We exclude it-band photometry due to its low 
S/N, but we note that this does not significantly effect 
the quality of our fits, as this data is of much lower qual- 
ity than in griz. 

In Figure [27] we present the SDSS-II light curves for 
the four SN candidates whose BOSS host-galaxy spec- 
tra are shown in Figure [31 We also present the SALT2 
best-fit SN la model and the one-sigma error on this fit 
provided by SALT2. 
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Fig. 27. — The r-band light-curves fitted using SALT2 for the SNe whose host galaxy spectra are shown in Figure [3] The blue solid 
curve is the best-fit SN la model light curve, and the red dashed lines represent the lcr uncertainties on this fit. The vertical dotted line 
shows the best— fit time to of peak brightness. 



C. REDSHIFT COMPARISON 

In Figure [28] we compare the BOSS host-galaxy 
rcdshift and the SN spectroscopic redshift for the 186 
spectroscopically-classified SNe la that pass our se- 
lection cuts. In only four cases do these redshift 
measurements disagree significantly (SN3199; SN13956; 
SN15301; SN19757). Based on visual inspection of the 
SDSS-III DR8 catalogue we identify and remove SN3199 
from our sample, as the BOSS-targeted galaxy does not 
appear to be the most likely SN host galaxy. The remain- 
ing three SNe are retained in our sample as the identified 
host galaxy appears to be correct and they reside close to 
the Hubble diagram (small HRs) when using the BOSS 
host-galaxy redshifts rather than the SN-spectrum red- 
shift. Apart from this check, the purpose of which is 
to determine the likelihood of host mis-identification, we 
make no other use of SN spectroscopy in this paper. 

We have also compared our BOSS host-galaxy red- 
shifts with SDSS galaxy redshifts, where the latter are 
available. We found one SNe la (SN6491) hosted in a 
galaxy where the redshifts from these two surveys dis- 
agree, and removed this object from our sample. 

Additionally, we have visually inspected the host 
galaxies of other photometrically-classified SNe la candi- 
dates in our sample. First, we inspected a random subset 
of 70 photometric SNe la, finding no obvious misiden- 
tification of the appropriate host. This is reassuring, 
as it confirms that the rate of misidentification of hosts 
must be low. Next, we inspected the host association for 
SNe la candidates that are clearly offset from the Hub- 
ble diagram in Figure [12] to ensure the correct galaxy had 
been assigned during targeting. We found only one host 
galaxy that was likely to be incorrect (SN9052), being lo- 
cated 25 arcseconds from our photometrically-classified 
SN, and removed this object from our sample. In total 



we removed only two SNe la from our sample because of 
likely host-galaxy mismatches. 

D. OTHER SYSTEMATIC UNCERTAINTIES 

There are a number of other systematic uncertain- 
ties that are likely to affect our SDSS-II SN sample be- 
yond the photometric-classification specific uncertainties 
discuss ed in Section 17.11 As ou tlined in iKessler et al.1 
(|2009al) and iConlev et all (|201lD . present SN samples 
have major uncertainties associated with their photomet- 
ric calibration and the light-curve fitting technique used, 
as well as many astrophysical uncertainties such as corre- 
lations with host-galaxy properties, SN lensing, peculiar 
velocities, galactic dust, and SN evolution. Though we 
do not address these important systcmatics in detail in 
this paper, we discuss the likely effect of these additional 
systematic uncertainties on our results. 

The optimal method of light-curve fitting is not 
known; there can be significant differences in the cos- 
mological results obtained from using different algo- 
rithms (e.g., SALT2 an d MLCS2 k2; IKessler et al.ll2009at 
Sollerman etaTl 120091) . Both IGuv et al.l (|2010f ) and 



Conlev et all (120111) find consistent cosmological results 



between SALT2 and a different technique (SiFTO), with 
a possible systematic uncertainty of only A// ~ 0.02 — 
0.03 magnitudes between light-curve fitters. We have 
chosen to use only the SALT2 light-curve fitting algo- 
rithm (|Guv et al.ll2010[) in this work . 

As discussed in [Conlev et all (|201lD for SNLS, 
the most important systematic uncertainty in present 
SN surveys i s the p hotometric calibration. Therefore, 
IConlev et al.l (|2011D recommended that future SN sur- 
veys should be calibrated onto a "more modern, better 
understood photometric systems such as USNO/SDSS'\ 
By using photometric data obtained wholly from the 
SDSS-II SN Survey, we believe we have minimized cal- 
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Fig. 28. — Comparison of the spectroscopic redshift for our spectroscopically-confirmed SNe la and the corresponding host-galaxy redshift 
from BOSS. The bottom panel shows the redshift residuals (BOSS galaxy redshift - SNe Redshift) for this sample. 



ibration uncertainties, as the SDSS phot ometric sys- 
tem i s now mature a nd well-understood (jlvezic et al.l 
l2007tlDoi et~aLll20T0h . iMosher et all (|20ll has recently 
compared the SDSS-II photometric system to that of 
the Carnegie Supernova Project (CSP) using light-curve 
data for nine SNe la observed concurrently by the two 
projects. They conclude that measurements from the 
two surveys agree in all bands at or better than 2% in 
flux, and are consistent with no difference in g— and 
r— band magnitude at the 2tr level. This is an indication 
of the relative calibration between the SDSS and other 
surveys (CSP in this case) and not a direct statement 
on the absolute calibration, which would require obser- 
vation of a known source like NIS T photodiodes, as dis- 
cussed bv lStubbs fc Tonrvl (|2012| ). Such techniques will 
be implemented in future surveys (e.g., DES). Based on 
these findings, we assume systematic uncertainty of only 
AyU ps 0.02 mags on the distance modulus of our SNe, 
relative to other surveys. 

The next largest systematic uncertainty at present 
is the recently discovered correlation between the (cor- 
rected) peak absolute magnitude of SNe la and the stellar 



mass of the host galaxy dKellv et al.H2010l : lSullivan et al.l 
l2010t iLampeitl et al.ll2010bl ). This relationship has been 
observed to have an effec t as large as Ap, ~ 0.07 mags in 
the SDSS-II SN sample (ILampeitl et alJl2010bl) . We do 
not attempt to make this correction here, as it is beyond 
the scope of this paper, but we note it is being inves- 
tigated in Smith et al. (in prep.). There will likely be 
significant stellar population modeling uncertainties as- 
sociated with determining reliable stellar masses of our 
faint host galaxies; the large magnitude errors on the 
photometric colors of our BOSS galaxies make this es- 
pecially difficult at z > 0.3. We will revisit the topic of 
SN la correlations with host-galaxy properties using our 
photometrically-classified sample in a future paper. 

There are a number of other astrophysical uncer- 
tainties that could be considered, especially SN lens- 
ing effects, peculiar velocities and possible SN evolution. 
These are all smaller in size compared to the system- 
atic uncertainties discussed above, and should be further 
mitigated in our sample because of the lack of relatively 
high-rcdshift SNe la where these effects are most promi- 
nent (especially lensing and evolution). 
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Fig. 29. — The 68% and 95% confidence contours of w and f2 m 
in a flat toCDM model. We assume a Gaussian prior on Hq; other 
priors are summarised in Table [5] using prior Set III. The blue 
contours show the statistics only contours and the grey contours 
include the estimate of the systematics, as well as statistical un- 
certainties. 



E. PHOTOMETRICALLY-CLASSIFIED SNE IA DATA 

In Table [6] we present the key information used 
in this paper for our sample of 752 photometrically- 
classified SNe la. The full SDSS-II SN sample, including 
all the light-curve data, rcdshifts and classifications for 
all transients, will be published in Sako ct al. (in prep.). 
Table [6] includes a unique identification number for the 
whole SDSS-II SN Survey (CID; column 1), BOSS host- 
galaxy rcdshift and error (columns 2 & 3), the RA and 
DEC (in degrees) of the SN event (columns 4 & 5) and 
its host (columns 6 & 7), a unique SDSS object identifier 
for the host galaxy from DR8 (column 8), the SALT2 
parameters Xq in flux units, X\ and color both in mag- 
nitudes (columns 9, 10, 11), and finally the uncorrected 
and corrected distance modulus (columns 12 & 13) with 
error (column 14), all in magnitudes. 

The data in Table H along with the SALT2 
covariancc matrices and the SN type proba- 
bilities, can be electronically downl oaded from 
http ://www. icg.port . ac .uk/^campbelh] The 
probabilities listed there are those used in this paper 
(with the BOSS host galaxy redshift prior) and include 
the probability assigned to each object of being a SNe la 
(P Ia ), a Type II SN (P 77 ), and a Type Ibc (P Ibc ). 



In Figure [29] we show an estimate of the effect of 
systematic uncertainties on our results, compared to Fig- 
ure [25] in Section 16.31 which includes only statistical er- 
rors. These measurements are obtained by adding in 
quadrature an additional uncertainty of A/i = 0.1 mag- 
nitudes to the distance moduli errors of our observe d 
SNe la; the same methodology as iKowalski et all (|2008h . 
This level of uncertainty is an estimate from the combi- 
nation (in quadrature) of uncertainties associated with 
the light-curve fitter (0.03 mags), photometric calibra- 
tion (0.02 mags), lensing and peculiar velocities (0.05 
mags) and possible host galaxy correlations (0.07 mags). 
When including our systematic errors in the cosmolog- 
ical analysis we fix the SALT2 SN parameters a and /? 
to the values found in Section [6] this is done to prevent 
a and j3 changing values to counteract the increase in 
dispersion that the systematic errors cause. The results 
do not change significantly with our inclusion of these 
estimates of the unknown systematics, giving a best-fit 
of w = — 0.93±g || and fl m = 0.314±g;g|. However, we 
stress again that this analysis is not comprehensive. 
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TABLE 6 

Table containing data for 10 of the 752 photometri cally-classified SNe Ia presented in this paper. For the full table, which in cludes errors and covariance on 

THE SALT2 PARAMETERS, SEE THE ELECTRONIC TABLE AT |hTTP://wWW.ICG.PORT.AC.UK/STABLE/CA MPBELH/SDSS-PhQTOME TRIC-SNe_IA.FITS| THE 15 SNE WITH ENTRIES OF * IN 
THE ELECTRONIC TABLE ARE ONES WHERE THERE IS NO PHOTOMETRIC OBJECT ID FOR THE HOST GALAXY IN DK8; THESE GALAXIES DO APPEAR IN THE CO- ADDED IMAGES, AND FROM 

THIS CATALOG WE QUOTE THE HOSTlD. 
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